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QUALITATIVE BEHAVIOUR OF INCOMPRESSIBLE 

TWO-PHASE FLOWS WITH PHASE TRANSITIONS: 

THE CASE OF NON-EQUAL DENSITIES 

JAN PRUSS, SENJO SHIMIZU, AND MATHIAS WILKE 



Abstract. Our study of a basic model for incompressible two-phase flows 
with phase transitions consistent with thermodynamics in the case of constant 
but non-equal densities of the phases, begun in I23| is continued. We extend 
our well-posedness result to general geometries, study the stability of the 
equilibria of the problem, and show that a solution which does not develop 
singularities exist globally. And if its limit set contains a stable equilibrium 
it converges to this equilibrium as time goes to infinity, in the natural state 
manifold for the problem in an Lp-setting. 
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1. Introduction 

Let ft C M" be a bounded domain of class C"^^ , n > 2. ft contains two phases: at 
time t, phase i occupies subdomain r2i(i) oiil. Assume dfli{t)nd^ = 0; this means 
no boundary intersection. The closed compact hypersurface T{t) := dCli(t) C fl 
forms the interface between the phases. 

Let u denote the velocity field, tt the pressure field, T{u, tt, 9) the stress tensor, 
D{u) = (Vm -|- [Vm]^)/2 the rate of strain tensor, 9 the (absolute) temperature 
field, vr the outer normal of f2i, Mr the velocity field of the interface, Vp = ur • i^r 
the normal velocity of r(t), Hr = H{r{t)) = — divri^r the curvature of r(i), j the 
phase fiux, and |w] = W2 — f i the jump of the variable v across T{t). 
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2 J. PRUSS, S. SHIMIZU, AND M. WILKE 

Several quantities are derived from the specific free energies ipi{d) as follows. 
ei{6) := ^i{9) + 9rji{9) means the internal energy in phase i, ?7i(0) = — V'K^) the 
entropy, Ki{e) := e[{6) > the heat capacity, and l{9) = e{iP'{e)\ = -9{ri{e)\ 
the latent heat. Further di{0) > denotes the coefficient of heat conduction in 
Fourier's law, iJ,i{d) > the viscosity in Newton's law, pi,p2 > the constant, 
positive densities of the phases, and a > the constant coefficient of surface 
tension. In the sequel we drop the index i, as there is no danger of confusion; we 
just keep in mind that the coefficients depend on the phases. 

By an Incompressible Two-Phase Flow with Phase Transition we mean the fol- 
lowing problem with sharp interface: 

Find a family of closed compact hypersurfaces {T{t)}t>o contained in il and ap- 
propriately smooth functions u : M.^ x fj ^- K.", and tt, 6* : IR+ x f2 ^- M such 
that 



p{dtu + M • Vm) - div T(u, TT,e) = 


in n\T{t), 




T{u,n,e)^2p{e)D{u)~nI, divu^O 


in n\T{t), 




-Jj'^r - lT{u, 7T, 9)i^TJ = (rHrur M = {-\3^t 


on r(i). 


(1.1) 


u = Q on 9ri, U|j^Q = uo 


inVL\T[t). 




'i{9){dte + u • V6I) - div {d{9)W9) - 2p\D{u)\l = 


in n\V{t), 




euow - ld{e)d,rO\ = o, M = o 


on r(t), 


(1.2) 


d^e = on dVL, 6i|,^„ = 9^ 


in VL. 




1^(^)1 --[^^.r + f^-'^'^^^^-^^l 


on r(i). 




Vt — UT ■ J^v — u- vr j 

p 


on r(i). 


(1.3) 


Tu.o - To. 







This model is explained in more detail in our previous paper |22| . It is based on a 
model recently proposed by Anderson, Gurtin et al. [2], see also the monographs 
by Ishii [11] and [12], and it is thermodynamically consistent in the sense that in 
absence of exterior forces and heat sources, the total energy is preserved and the 
total entropy is nondecreasing. It is in some sense the simplest sharp interface 
model for incompressible Newtonian two-phase flows taking into account phase 
transitions driven by temperature, which is consistent with thermodynamics. We 
want be more precise on this topic. 

The total energy of the system is given by 



£{u,d,T)^ [ [t\u\^ + pe]dx + cj\Vl 
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and the total entropy reads 



$(6»,r) = / pr]{e)dx. 



For the case pi ^ p2 the following result is valid. 

Theorem 1.1. Let a^p\,p2 > 0, pi 7^ p^, and suppose that tpj € C^(0,oo), 
Pj,dj e C'^{0,oo), 

-i;]is),pj{s),dj{s) > 0, for s > 0, j = 1,2. 

Then the following assertions hold. 

(i) The total energy E is conserved along smooth solutions. 

(ii) The negative total entropy —^ is a strict Ljapunov functional, which means 
that —^ is strictly decreasing along nonconstant smooth solutions. 
(iii) The equilibria without boundary contact are zero velocity, constant tempera- 
ture, constant pressure in each phase, and f^i consists of a finite number of non- 
intersecting balls of equal size. 

(iv) The critical points of the entropy functional with prescribed energy and phase 
volumes are precisely the equilibria of the system. 

(v) Local maxima of the entropy functional with prescribed energy and phase vol- 
umes are precisely the equilibria with connected interface. 

In equilibrium, the size of the balls follows from eonservation of total mass 

Pl|r2i| + P2\^2\ = Co, 

hence 

mK/n)i?" = (p2|^|-co)/Ip]; 
once the number of components m of the disperse phase Jli is given, this determines 
the radius i? > of the balls. The temperature is uniquely determined by the total 
energy 

pi|r!i|ei(0) + P2\n2\e2{e) + (Tmw„i?"-i = Eo, 
since the functions ej{9) are strictly increasing by assumption. Therefore we have 
an n?7i-paramcter family of equilibria, where the parameters are the centers of the 
TO balls. 

It turns out that in case m = 1 each equilibrium is linearly stable and as 
one main result in this paper we prove its nonlinear stability. We also prove 
on the other hand, if m > 2 then each equilibrium is unstable, linearly as well 
as nonlincarly; this is a well-known phenomenon called Ostwald ripening in the 
literature. Thus local maxima of the entropy under the constraints of given energy 
and phase volumes correspond to the stable equilibria of the system. This is 
precisely what one would expect physically. 

There is a large literature on isothermal incompressible Newtonian two-phase 
flows without phase transitions, and also on the two-phase Stefan problem with 
surface tension modeling temperature driven phase transitions. On the other hand, 
mathematical work on two-phase flow problems including phase transitions are 
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rare. In this direction, we only know the papers by Hoffmann and Starovoitov 
PITO] dealing with a simplified two-phase flow model, and Kusaka and Tani p^ITT] 
which are two-phase for temperature but only one phase is moving. The papers of 
Di Benedetto and Friedman [5^ and Di Benedetto and O'Learyifij consider weak 
solutions of conduction-convection problems with phase change. However, none of 
these papers deals with models which are consistent with thermodynamics. 

We emphasize that the major difference between the cases of equal or different 
densities lies in the occurrence of the so-called Stefan currents which arc induced 
by the jump in the normal velocity across the interface in case pi y^ P2- If the 
densities are equal these are absent, this case which we call temperature- dominated 
has been analyzed in [22j and |27j . Here we are interested in the velocity-dominated 
case, i.e. densities pi ^ P2. Note that in this case we may eliminate the phase flux 
j, multiplying the jump condition for u by i^r, to the result 

j = |,i.i/r]/|l/pl. 

The jump condition for u then becomes -PtIm] = 0, where Pr = I — vr 'Si vr- 
Similarly, multiplying the second equation in ()1.3|) with p and taking the jump 
across T{t) we obtain the new equation 

IpjVr - Ipu ■ z.r] 

for the moving interface r(i). However, it is useful at several occasions to keep 
the variable j, as it has a clear physical meaning. 

In [23] the first two authors investigated local well-posedness of this problem 
in a situation of a nearly flat interface, based on Lp-maximal regularity of the 
underlying principal linearization. It is the purpose of this paper to extend this 
work to general geometries and to present a qualitative analysis of (jl.ip , (|1.2p , (|1.3|) 
in the framework of Lp-theory. We study the equilibria of the system which are zero 
velocities, constant pressures in the phases, constant temperature, and the disperse 
phase Qi consists of a finite number of non-intersecting balls of the same radius. 
Thus the equilibria form a manifold, and we prove that an equilibrium is normally 
stable if and only if the phases are connected, otherwise it is normally hyperbolic; 
see below for a definition of these notions. We further prove by means of the 
generalized principle of linearized stability that an equilibrium is stable for the 
nonlinear problem if and only if the phases are connected. Furthermore we show 
that a solution which does not develop singularities, in a sense specified below, 
exist globally in the natural state manifold SA4 of the system, and it converges 
to a probably different equilibrium, provided its limit set in SM contains a stable 
equilibrium. 

The plan for this paper is as follows. In Section 2 we employ the direct mapping 
approach to transform the problem locally in time to a fixed domain. This is by 
now a well established method, we refer in particular to [24j for the necessary geo- 
metric background. Section 3 deals with local well-posedness, which is based on 
maximal Lp-regularity of the underlying principal linearization and the contraction 
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mapping principle. The proof of the crucial maximal regularity result for the aris- 
ing non-standard asymmetric Stokes-problem is given in Section 7; by perturbation 
and localization, it is derived from the corresponding result for flat interfaces [2 3) . 
In Section 4 we study the linearization of the problem at a non-degenerate equi- 
librium. Employing the generalized principle of linearized stability, in Section 5 
we prove these stability assertions also for the nonlinear problem. In Section 6 we 
introduce the natural state manifold SA4 of the system in the Lp-setting, and it is 
shown that the local well-posedness result from Section 3 induces a local semiflow 
on SA4. We study its asymptotic behavior, employing the negative entropy as 
a strict Ljapunov functional, relative compactness of bounded semiorbits and the 
stability result from Section 5. 

In a subsequent paper we will extend the results of this contribution to the 
case where the coefficient of surface tension ct > is no longer constant but 
temperature-dependent. This allows us to include the so-called Marangoni forces. 

2. Transformation to a Fixed Domain 

Let Q, C K" be a bounded domain with boundary dO, of class C^ , and suppose 
r C fi is a hyper-surface of class C^, i.e. a C^-manifold which is the boundary 
of a bounded domain Jli C lli C fi; we then set 0,2 = ri\f2i. Note that VI2 
is connected, but Vii maybe disconnected, however, it consists of finitely many 
components only, since dQ.i = F by assumption is a manifold, at least of class C^. 
Recall that the second order bundle of F is given by 

N^V := {(p, vv{p),Vtvt{p)) : P e F}. 

Here vt{p) denotes the outer normal of J7i at p e F and Vr the surface gradient 
on F. The Weingarten tensor Lt on F is defined by 

Lt{p) = -Vrvrip), p G F. 

The eigenvalues Kj{p) of Ly{p) are the principal curvatures of F at p e F, and we 
have |Lr(p)| = niaxj |kj(p)|. The curvature Hr(j)) (more precisely (n — 1) times 
mean curvature) is defined as the trace of Lt{p), i.e. 

n-l 

■ffr(p) = ^Kj ^ -divriyrip), 

where divr means surface divergence. Recall also the Haussdorff distance dn 
between the two closed subsets A,Bc M™, defined by 

dniA, B) := maxjsup dist(a, B), sup dist(6. A)}. 

Then we may approximate F by a real analytic hyper-surface E (or merely S G C^), 
in the sense that the Haussdorff distance of the second order bundles of F and E is 
as small as we want. More precisely, for each 77 > there is a real analytic closed 
hyper-surface E such that d//(A/'^E,A/'^F) < 77. If 77 > is small enough, then E 
bounds a domain ilf with Vif C fi, and we set J7f = il\rjf . 
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It is well known that a hyper-surface S of class C^ admits a tubular neighbor- 
hood, which means that there is a > such that the map 

A : S X (-a, a) ^ R" 
A{p,r) := p + rvj:{p) 

is a diffeomorphism from E x {—a, a) onto TZ{A), where I'sip) denotes the outer 
normal at p € S. The inverse 

A"^ :7^(A) K^ E X {-a,a) 

of this map is conveniently decomposed as 

Here Il-s{x) means the orthogonal projection of a; to E and d-s{x) the signed 
distance from x to E; so |(is(x)| — dist(x, E) and ds{x) < if and only if a; € ilf. 
In particular we have Tl{A) — {x E R" : dist(x, E) < a}. 

Note that on the one hand, a is determined by the curvatures of E, i.e. we must 
have 

< a < min{l/|Kj(p)| : j = I, . . . ,n — 1, p £ E}, 

where Kj{p) mean the principal curvatures of E at p G E. But on the other hand, 
a is also connected to the topology of E, which can be expressed as follows. Since 
E is a compact manifold of dimension n — 1 it satisfies the ball condition, which 
means that there is a radius r^ > such that for each point p G E there are 
Xj S flf, j — 1,2, such that Br^{xj) C rij', and Br^{xj) n E = {p}. Choosing 
rs maximal, we then must also have a < r-^:. Note that l/rs also bounds the 
principal curvatures, we have 

kjb)| < l/^E, j = l,...,n- 1, pe E. 

In the sequel we fix a = r-s/2. 

For later use we note that the derivatives of lis (x) and ds (x) are given by 

Vds(x) = ^s(ns(a:)), Z?(u)ns(x) = A/o(ds(.T))(ns(a;))Ps(ns(x)), 

where Py,{p) — I — vy,{p) ® v^{p) denotes the orthogonal projection onto the 
tangent space T^E of E at p e E, and Mo(r)(p) = (/ — ris(p))~^- Note that 

|Mo(r)(p)| < 1/(1 - r\L^{p)\) < 2, for ah p e E, \r\ < a. 

All of these facts are discussed in more detail in [24] . 

Setting r = r(i), we may use the map A to parameterize the unknown free 
boundary T(t) over E by means of a height function h{t,p) via 

r{t):p^p + h{t,p)uj:{p), peE, 

for small i > 0, at least. Extend this diffeomorphism to all of 0, by means of 

Eh{t,x) = X + xid^{x)/a)h{t,Il^{x))v^{Il^{x)) =: x + ^h{t,x). 
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Here x denotes a suitable cut-off function; more precisely, x € 2?(]R), < x < 1, 
X{r) = 1 for \r\ < 1/3, and x(^) = for \r\ > 2/3. Note that Eh{t,x) = x for 
Me (2;) I > 2a/3, and 

'^h^{t,p)=p-h{t,p)vT.{v) for pes. 

Now wc define the transformed quantities 

u{t,x) = u{t,Eh{t,x)), 9{t,x) = d{t,Ehit,x)), 

j{t,x) =j{t,Eh{t,x)), Tf{t,x) ='K{t,Eh{t,x)), t>0,xefl\T,, 

the pull backs of u, 9, j, and tt. This way we have transformed the time varying 
regions fi \ r(t) to the fixed domain fi \ E. 

This transformation gives the following problem for u, tt, 9,j, h. 

pdtu - g{h) ■ i^{9)(g{h)u + [g{h)u]'^) + g{h)7T = 7^„(^2, 9, h) in f7\S, 

g{h) • u = in r2\s, 

pK{9)dt9 -g{h)- d{9)g{h)9 = Tie (u, 9, h) in f7\S, 

u — d^9 — on 9ri, 

[l/pffiyrih) + l~^li9)igih)u + [g{h)u]^) + Ttjiyrih) = aHr{h)iyrih) on S, 

PrN^O, J=[w-H/[1/Pl, 1^1=0 onS, (2.1) 

Ii9)j + ldi9)gih)9-iyr{h)]=0, on S, 

1^(0)1 + Il/2p2]P _ |_(^(0-)/p)(g(;i)s -^ [g{h)uV) - 7T/pPr{h) = on S, 

Ipl^r - Ipu ■ i^rj = 0, on S, 

u(0) = uo, e(0) = 00, h{0) = ho. 

Here 5(ft.) and Hr{h) denote the transformed gradient and curvature, respectively. 
More precisely, we have 

DEh^I + D^h, [DEh]-' =!-[! + DihV^Dih =1- M^(hY . 

A simple computation yields 

D^h{t,x) = u^{Jl^{x)) ® Mo{d^{x)){Jl^{x))W^h{t,W^{x)) 

-/i(i,ns(a;))Ls(Hs(a;))Mo(rfs(a;))(ns(.T))Ps(ns(a;)) 

for |rfs(2;)| < a/3, DS,h{t,x) ~ for |(is(a;)| > 2a/3, as well as 

D£,h{t,x) = -x'{d^{x)/a)h{t,IV^{x))v^{I{^{x)) ® MlV^{x)) 
a 

+ x(ds(a:)/a)[i.E(HE(x))®Mo(ds(a:))(HE(x))Vs/i(i,Hs(x)) 

-/i(i,Hs(x))Ls(Hs(x))Mo(ds(x))(Hs(x))PE(Hs(x))] 

for a/3 < Ms(x)| < 2a/3, 

Therefore [/ -I- D£^h] is invertible, provided h and V^ft. are sufficiently small, more 
precisely 

\[I + D^nV\ < 2, for \h\^ < 1 min{a/|x'U, VlisU}, \V^h\^ < ^. (2.2) 
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With these properties we get 
[Vtt] o Eh = g{h)7t 

= {[DEh]-^ o Eh)Vn - Vtt - [D^hV[I + D^hr^Vit 
=: {I-Mi{h))VTt 

[ve]oEh = {i-Mi{h))ve 

divuoEh = g{h) -u = (I - Mi{h))V ■ u. 
Next we note that 

dtu oEh^ dtu - [Vu] o Eh ■ dtEh = dtu - {[DEh]~'^ o Eh)Vu ■ dtEh 
= dtu -Vu-[I + DChr^dt^h =■■ dtu - (Rih) ■ V)u, 
with R{h) ^{I - Mi{hy)dtS,h- This yields 

7^„(u, e, h) = -pu ■ g{h)u + p{R{h) ■ V)u, 
and 

ne{u, 0, h) = -pK{e)u ■ g{h)0 + pK{e){R{h) ■ v)0 
+ fiie){g{h)u + [g{h)u]'^) ■ g{h)u. 

With the Weingarten tensor Ls and the surface gradient Vs we further have 
vr{h) = I3{h){u^ - a{h)), a{h) = Moih)Vj:h, 

Mo{h) - (/ - hL^)-\ P{h) = (1 + \a{h)\Y'^\ 

and 

Vr = {dtE\iyr) = dth{,^rW^) = f3{h)dth. 
The curvature Hr{h) becomes 
Hrih) = /3{h){tT[Moih)iLj: + Vj:aih))] - l3^h)iMoih)aih)\[V^aih)]a{h))}, 

a differential expression involving second order derivatives of h only linearly. Its 
linearization at /i = is given by 

i/f(0)=trL| + As. 

Here As denotes the Laplace-Beltrami operator on S. 

It is convenient to decompose the stress boundary condition into tangential and 
normal parts. Multiplying the stress interface condition with i^e//? we obtain 

[1/pl.f + 1^1 - 'yHrih) - {lp{9){g{h)u + [g(/i)S]T)](j.s - Mo{h)^sh)Ws), 
for the normal part of the stress boundary condition, and 

- PJ:l^l{e){g{h)u + [g{h)u]'')j{i^^ - Mo(h)w^h) 

= (im(g{h)u+[g{h)u]'')j{iy^-Mo(h)W^h)\iy^)Mo{h)W^h, 



for the tangential part. Note that the latter neither contains the pressure jump 
nor the phase flux nor the curvature, which is the advantage of this decomposition. 
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3. Local Well-Posedness 

In this section we prove local well-posedness of problem (|l.l|) . (|1.2|) . (|1.3p . The 
proof is based on maximal Lp-regularity of the principal linearization and on the 
contraction mapping principle. 

3.1. Principal Linearization. The principal part of the linearized problem reads 
as follows 

pdtu — fi{x)Au + Vtt = pfu in il \ E, 

div u = gd in il \ S, 

PsN + c(t,a;)Vs/i = ftSn on S, (3.1) 

-2ln{x) D{u)v^\ + {-kP^ - o-A^hvj: =5 on S, 

u = on dfl 

w(0) = uq in rj \ S. 

pK{x)dte - d{x)A0 = pK{x)fe in f7 \ E, 
-[d{x)dv^e} ^ ge on S, 

|6i] = on S, (3.2) 

d^e = on 9f7 
61(0) ^ 00 in n. 

Ipjdth - Ipu ■ J^sl + Kt, x) ■ V^h = Ipjfh on S, 
-2lif,ix) /p)Diu)i^^ ■ i/sl + [tt/p] = gh on S, (3.3) 

/i(0) = ho on E. 

Here pk,K-k,dk, k = 1,2, are functions of x, continuous on fjfc, and c, 6 depend 
on t and x; recall that |p] 7^ by assumption. Apparently, p.2p decouples from 
the remaining problem. Since it is well-known that this problem has maximal 
Lp-regularity, we concentrate on the remaining one for (M,7r, /i). 

pdtU — p,{x)Au + Vtt = pfu in i^ \ S, 

div u = gd in ri \ S, 

u — Q on dVt 

PsN + c{t, x)Vsh = P^gu on E, (3.4) 

-21p{x)D{u)u^\ + IttJi^s - a-A^hiy^ =g on E, 

-2|^(x)D(u)i/s • i's/pI + [tt/p] = 5h on E, 

|pl9t/i - |pu • i^sj + b{t, x) ■ Ve/j = Mfh on E, 

m(0) = uo in rj \ E, 

/i(0) = ho on E. 

For this problem we have maximal regularity result in the Lp-setting. 
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Theorem 3.1. Let p > n + 2 be fixed, pj > 0, p2 ^ pi, fij G C{flj), pj > 0, 
J = 1,2, (&,c)€ [W^;-^/'P(J;Lp(E))nLp(J;WV'-^/^(I]))]"+i, where J ^ [0 , a] . 

Then the Stokes problem with free boundary p.4p admits a unique solution 
(u, TT, h) with regularity 

u e Hl{j- Lp{n)r n Lp{j; Hl{n \ S))", 

z/ and only if the data fu, gd, 9, P^gu, gh, fh, uq, Hq satisfy the following regularity 
and compatibility conditions: 

(a) /„eLp(J;Lp(r!,M"+i)), 

(b) gd e H^{J;H-\n)) n Lp{J;H^{n \ S)), 

(c) (g,g/.) eT</'-'/'^(J;Lp(E,M"+i))nLp(J;T<-^/^(I],M"+i)), 

(d) (P^gujh) G t^^'/'^(J;Lp(S,M"+i)) nLp(J;iyp'-^/P(E,M"+i)), 

(e) uo e iyp'-'/''(n \ E,K"), /lo e WV'"'/^(S), 

(f) divMQ = gd(0) m il \ S, 

(g) Ps W + c(0, OVs/io = ^s5«(0) on S, 

(h) -Pe[mo(-)(V«o + [Vuo]^)1=Ps5(0) on E. 
T7ie solution map [{fu,gdi <?, -Psffu, <?/i, A, wq, /iq) >-> (w, tt, /i)] is continuous between 
the corresponding spaces. 

The proof of this result is given in the Section 7. 

3.2. Local Existence. The basic resuh for local well-posedness of Problem (|l.ip . 
(|1.2p . (|1.3p in an Lp-setting is the following theorem. 



Theorem 3.2. Letp > n+2, fi C M" a bounded domain with boundary dfl G C^ , 
a, pi, P2 > 0, pi 7^ p2, and suppose 4'j G ^^"^(0, oo), pj, dj G C^(0, oo) are swc/i that 

Kj{s) ^ -sTpjis) > 0, pj{s)>0, dj{s)>0, s G (0, oo), j = 1,2. 

Assume the regularity conditions 

iuo,eo) G [w^~'/p{n\To)r+\ To G %3-2/p, 

the compatibility conditions 

divuo = in fl\ Fq, uq = O^Oq = on dfl, 
^roNl = Pro[M(^o)(V^io + [Vuo]^)^ro] - I^o] - on Tq, 
Z(0o)N • ^rol/Il/p] + Id(0o)5.r„"o] =0 on Tq, 
and t/ie well-posedness condition ^q > on Cl. 
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Then there exists a unique Lp-solution of Problem il.l\) . il.2\) . il.3\) on some 
possibly small but nontrivial time interval J = [0, r] . 

Here the notation Fq S Wp means that Fq is a C^-manifold, such that its 

(outer) normal field vyo is of class Wp~ (Fp). Therefore the curvature tensor 
Lr„ = -Vroi^Fo of Fo belongs to Wp~^^''{To) which embeds into C"+i/p(Fo), with 
a = 1 — (n + 2)/p > since p > n + 2hy assumption. For the same reason we also 
have Mo G Ci+"(f^j(0)))", and 6*0 e Ci+"(f^j(0))), j = 1,2, and Vq G C1+"(Fo). 
The notion Lp-solution means that (u, w, 9, F) is obtained as the push-forward 
of an Lp-solution {u,Tr,0,h) of the transformed problem (|2.ip . which means that 
{u, 9, h) belongs to E(J) = E„^0(J) x E;j(J) defined by 

¥.uAJ) = {(", 9) G [Hl{J:. Lp{n)) n Lp{J; H^{n \ E))]"+i : divu = 0}, 
and 

E,(J) := Typ2-1/2P(J;L^(S)) ni/pi(J;M/p2-i/P(E)) nLp(J;M^p3-i/P(E)). 
The regularity of the pressure is obtained from the equations. 

3.3. Time- Weights. For later use we need an extension of the local existence 
results to spaces with time weights. For this purpose, given a UMD-Banach space 
Y and /x e (1/p, 1], we define for J — (0, to) 

k;JJ; y) :- {u e LpjociJ; Y) -. t^-^u e K;iJ; r)}, 

where s > and K e iJ, VF. It has been shown in [T^] that the operator d/dt in 
Lp,fi{J',Y) with domain 

i?(d/dt) - oHiAJiY) ^{ue Hl^{J-Y) : u(0) = 0} 

is sectorial and admits an _ff°°-calculus with angle 7r/2. This is the main tool to 
extend Theorem 13.21 to the time weighted setting, where the solution space E( J) 
is replaced by 

E^(J)=E^,„,0XE^.;,(J), 

with 

E^,„,e(J) = {(w,0)e [i/pV(J;Lp(f^))nLp,^(J;Lfp2(f^\S))]"+i : divu = 0}, 
and 

E^,,(J) := wl^'/'^iJ- Lp{j:)) n Hl^iJ; w^'^/^i:)) n Lp^; w^''/^{i:)). 

The trace spaces for u, 9 and h for p > 3 are then given by 

iuo,9o) € [W^^-'/P{n \ E)]"+\ /lo G VFp2+'^-2/p(5.)^ 
/ii := 9t/i|,^„ e VFp2^-3/p(j.)^ (3_g) 

where for the last trace we need in addition /i > 3/2p. Note that the embeddings 

E^,n,e( J) -^ C{J; C\n,)r+\ Ep,„( J) -^ C{J; ^^+"(1])) n Ci(J; C^CE)) 
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with a = 1/2 — n/p > require /i > 1/2 + (n + 2)/2p, which is feasible since 
p > n + 2 hy assumption. This restriction is needed for the estimation of the 
nonhnearities. 

For these time weighted spaces we have the following result. 

Corollary 3.3. Letp> n + 2, fie (1/2 + (n + 2)/2p, I], <7,pi,p2 > 0, pi ^ p2, 
and suppose ip S C'^(0, oo), p,d E C^(0, oo) are such that 

Kj{s) ^ -sipj{s) > 0, pj{s)>0, dj{s)>0, s e (0, oo), j = 1,2. 

Assume the regularity conditions 

K, Oo) e [w^^-'^n^ \ ro)]"+\ To e w^+^^-^^p, 

the compatibility conditions 

divuo = m il \ E, uq — d,j9o ~ on d^, 

Pro Nl = ^ro b(eo)(Vwo + [Vuo]^)^ro] - I^o] - on To, 
mnuo ■ ^rol/Il/p] + ldi9o)d,r,uo] - on To, 

as well as the well-posedness condition 9q > on i}. 

Then the transformed problem (|2.ip admits a unique solution z = (u, 9, h) € 
E^(0, t) for some nontrivial time interval J ~ [0, r]. The solution depends contin- 
uously on the data. For each 6 > the solution belongs to E((5, r), i.e. regularizes 
instantly. 



3.4. Proofs. The proof of Theorem 13. 21 follows similar ideas as in the papers [14] 
for the two-phase Navier-Stokes problem without phase transition, [22] for the 
problem with phase transitions and equal densities and |23| for the case of non- 
equal densities and nearly flat interface. Therefore we refrain from presenting all 
details but instead concentrate on the main ideas; note that the spaces are the 
same as in [53j. 

The transformed problem is rewritten in quasilinear form, dropping the bars 
and collecting its principal linear part on the left hand side. We set po{x) = 
p{9q{x)), Koix) = K(6'o(a;)), doix) = d{9Q{x)), c{t,x) = -e^^*|uo • J^sl, and 
b{t,x) — e ^*|puo]. Then, eliminating j, the problems reads as follows 

pdtu — pq{x)Au + VTr = Fu{u,TT, 9, h) in rj \ S, 

dYvu = Gd{u,h) in f^ \ S, 

pKo{x)dt9-do{x)A9 ^ Fe{u,9,h) on f7 \ S, (3.7) 

u — du9 — Q on do., 

u(0) = MO, 9{Q) = 6io in f7 \ S, 



on 


x; 


on 


E 


on 


Y] 


on 


E 


on 


Y] 


on 


x; 


on 


E 


on 


E 



(3.8) 
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Pj:Iu] + cit,x)Vj:h = Gjiu,9,h) 
-Psbo(x)(V«+ [Vm]T),.s] = G'r{u,e,h) on 
^2l^l^{x)Vu}v^ ■ v^ + W - G^y)i = GT\u, e, h) on 

{p\dth - {pu ■ v^\ + b{t, x)-W^h = Fh{u, h) 

KQ) = ho 

The nonlinearities are given by 

Fu{u, 9, TT, h) = {p{9) - p{9a))Au + Mi{h)V7r - p{u ■ {I - Mi{h)) - R{h)-)Vu 
+ /x'(^)((/ - Mi{h))V9 ■ {{I - Mi{h))Vu + [{I - Miih))Vu]^)) 

- ^i{e){M2{h) : V^)u - p{9){Nh{h) ■ V)m + ^i{e)MA{h) : Vu, 
Gd{u,h) = Mi(/i) : Vu, 

Fg{u, 9, h) = p(k(0o) - >i{e))dt9 - {d{9) - d{9o))A9 - d{9)M2{h) : V^9 

+ d'{9)\{I - Mi{h))W9\^ - d{9)M3{h) ■ W9 - pK{9){R{h) ■ W)9 

- Pk{9)u- {I- Mi{h))V9 
+ p{e){{I - Nh{h))Vu + [{I - A'h{h))\7uf) : {I - Nh{h))Vu, 

G'r{u,9,h)^P^l{^l{9)-^l{9M^u+[^uf)M 

- Pj:Ip{9){Vu + [VuY)Mo{h)V^h} 

- PY.lp{9){Mi{h)Wu + [Mi{h)WuY){uj: - Mo{h)\/^h)] 
+ H9){{I - K'h)Vu + [{I - Mi)Vvf) 

Gl°-{u, 9, h) = l{p{9) - p{9o)){\/u + [\/uf)u,: ■ ^^s] 

- M9)iVu + [Vw]T)Mo(/i)Vs/^ • ly^j 

- liii9)iMi{h)Vu + [Miih)Vu]''){iy^ - Mo{h)V^h) ■ i^^ 
+ a{Hr{h)-A^h)-lu-uryil/pl 

Ge{9, h) = l{9)j - ldo{x)d,9} + ld{9){I - lVh{h))V9 ■ ^^rl, 
GUu,9,h)^-liji9)]-ll/2p%f 

+ 2l{fi{x)/p)d,uj - l{p{9)/p){Mi{h)\7u+[M,{h)Wu]'')i,r ■ i^rj, 
Fh{u,h) = {b{t,x) - lpMo{h)ul)Vj:h, 

j^lu- ^s]//3(/i)[l/p], i^r = /3(/i)(^E - AUh)V^h). 
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Here we employed the abbreviations 

Kh{h)^[D^hV[I + Dih]-^V^, 

M2{h) = Mi{h) + Mj{h) - Mi{h)Mj{h), 

Msih) = {I-Mi{h))dWM2{h), 

M^{h) = {{I - Ah{h))V)Ah{h) - [((J- Mi(/i))V)Mi(/i)]T. 

We prove local well-posedness of (|3.7p , (|3.8I) by means of maximal Lp-regularity 
of the linear problem (Theorem 13. II) and the contraction mapping principle. The 
right hand side of problem (|3.7p . (|3.8p consist of either lower order terms, or terms 
of the same order as those appearing on the left hand side but carry factors which 
can be made small by construction. Indeed, we have smallness of /iqi Vs/iq and 
even of V|^/io uniformly on E, because Fq is approximated by S in the second 
order bundle. 6 appears nonlincarly in i/;, k, /i, d, but only to order zero; hence e.g. 
the difference {fi{9{t)) — /i(6'o)) will be uniformly small for small times. 

We introduce appropriate function spaces. Let J — [0,a\. The solution spaces 
are defined by 

Ei(a) -.^ {u e H^{J; Lpinyy n Lp{J; H^{n\j:)y' : u = Oondn,}, 

E2ia):=Lp{J:H^{Q\j:)), 

Eiia) ■= {9 e H^{J; Lp{n)) n Lp{J; H^{n \ S)) ■.0^9 = on dn, {9} = 0}, 

E,{a) := W^-'/^P{J;Lp{j:))nH'p{J;W^-'/P{j:))nLp{J;W^''/P{j:)). 

We abbreviate 

E(a) := {{u,TT,TTj,9,h) G Ei(a) x £2(0) x £3(0) x £4(0) x £5(0)}, 

and equip Ej(a) (j — 1, . . . , 5) with their natural norms, which turns E(a) into a 
Banach space. A left subscript always means that the time trace at t = of the 
function in question is zero whenever it exists. 
The data spaces are defined by 



Fi(a 
F2(a 
F3(a 
F4(a 
F5(a 
¥e{a 
¥r{a 
Ma 



■.^Lp{J;Lp{n)r, 

:= H;,{J; Hp\n)) n Lp{J; H;,{n \ E)), 

:- Ps[W^''/^P{J;Lp{j:;M.^))nLp{J;W^-'^P{^,R"))] 

:= W-V2-i/2p( J. ^^(5^))n n Lp{J; W^-^/Pi^)r, 

■.= LpiJ;Lpin)), 

:= W-V2-i/2p( J. ^^(2)) n LpiJ; W^~'^P{^)), 

:= W;i/2-i/2p( J. ip(s)) n Lp( J; W-pi-i/''(E)), 

:= W^,i-i/^P( J; ip(E)) n L,( J; ^^-^/^(E)). 
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We abbreviate 

8 

^(a) ■■= {{fu,gd,9j,9u,fB,9e,9h,fh) e ]jFfe(a)}, 

fe=i 

and equip Ffe(a) (fc = 1, . . . , 8) with their natural norms, which turns F(a) into a 
Banach space. 

Step 1. In order to economize our notation, we set z — (u, tt, iti,tt2, 0, h) G E(a) 
and reformulate the quasilinear problem p.7|) . p.Sp as 

Lz^N{z) (M(O),0(O),/i(O)) = K,0o,/io), (3.9) 

where L denotes the linear operator on the left hand side of p.7p . p.Sp . and 
N denotes the nonlinear mapping on the right-hand side of (|3.7I) . (|3.8p . From 
Theorem 13.11 we know that L : E(a) — ^ F(a) is bounded and linear, and that 
L : oF(a) — > oF(a) is an isomorphism for each a > 0, with norm independent of 
< a < flo < oo. 

Concerning the nonlinearity N ^ we have the following result. 

Proposition 3.4. Suppose p > n + 2, a > 0, and let ijji £ C^{0,oo), fii,di G 
C2(0, oo), K,{s) = -sip"{s), for i = 1, 2. 

Then for each a > the nonlinearity satisfies N € C"'^(E(a), F(a)) and its Frechet 
derivative N' satisfies in addition iV'(u, tt, tti, 7r2, 0, /i) € y8(oE(a), oF(a)). More- 
over, there is rj > such that for a given 2* G E(ao) with |/io|c2(E) !i V) there are 
continuous functions a{r) > and l3{a) > with a(0) = /3(0) = 0, such that 

\\N'{z + z*)||e(oE(a),oF(a)) < ckW + /3(a), z e Br C oF(a). 

This proposition is proved estimating the nonlinearities in the same way as in 
[23 and dg. 

Step 2. We reduce the problem to initial values and resolve the compatibilities 
as follows. Thanks to Proposition 4.1 in [T3], we find extensions gj g F2(a), 
5* S F3(a) which satisfy 

3^(0) = divwo, ft5:(0) = -Ps[/io(V7io + [VuoVy^l 

Further we define 

9* :=e^-*G,(0o,/^o), 9*9 := e^-'Ge{uo,0o,ho), 

and set f* ^ fg = ft = 9*h — ^- With these extensions, by Theorem 13.11 we 
may solve the linear problem (|3.1I) . (|3.2I) . (|3.3I) with initial data (uo,^Oi^o) and 
inhomogeneities (/* , gj , g* , 9u, f'e ^ 9*e , 9hJh) ^ '^hich satisfy the required regularity 
conditions and, by construction, the compatibility conditions, to obtain a unique 
solution 

z* = (w*,7r*,7r*,7r;,6l*,/i*) G E(a) 
with u*(0) = uo, 6'*(0) = 0Q, and /i*(0) = /iq- 
Step 3. We rewrite problem (13. 9p as 

Lz^ N{z + z*)- Lz* ^-.Kiz), zeoF(a). 
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The solution is given by the fixed point problem z = L^^K{z), since Theorem 13. II 
implies that L : oIE(a) — > o]F(a) is an isomorphism with 

|i~^l£(oF(Q),oE(a)) < M, a e (0,ao], 

where M is independent of a < gq. We may assume that M > 1. Thanks to 
Proposition 13.41 and due to K{0) — N{z*) — Lz*, we may choose a e (0, aq] and 
r > sufficiently small such that 

r 1 

|^(0)|F(a) < ^, |^'(z)|£(oE(a),oF(a)) < ^, Z G oE(a), \z\E{a) < r 



hence 



|i^WlF(a) < ^, 



which ensures that L~^K{z) : M^E(^a)iO,f) — > B(jE(a)(0,r) is a contraction. Thus 
we may employ the contraction mapping principle to obtain a unique solution on 
a time interval [0, a], which completes the proof of Theorem 13.21 CoroUarv 13.31 is 
obtained in the same way. 



4. Linear Stability of Equilibria 

1. We call an equilibrium non- degenerate if the balls do neither touch each 
other nor the outer boundary; this set is denoted by £. To derive the full lineariza- 
tion at a non-degenerate equilibrium e* := (0,0*,!]) £ £, note that the quadratic 
terms u ■ Vu, u ■ V6', |£'(u)||, |u]j, and |l/2p^]j^ give no contribution to the 
linearization. Therefore we obtain the following fully linearized problem. 



pdtu — fj,^,Au + Vtt ~ pfu in 57 \ E, 



where /i* = /i( 
temperature d 







div u^ ga 


in 
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s, 
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on 
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with K* = k{0*), d* — d{0*) and /* = 1(0*). The remaining conditions on the 
equilibrium interface S are 

|-p"^T(u, TT, i?)z^s • z^s] + Z*i9 = gh on S, 

9t/i-|pu-^sl/M-A onE, (4.3) 

h{0) = ho on E. 

The time-trace space E-, of E( J) is given by 

K, t?o, /lo) e E^ = [%'-'/''((f^ \ Er+^ X Typ3-2/P(E), 
and the space of data is 

{{fuJe),gd,{fh,PT.9u),{9,9e,9h)) e ¥{J) 

:= F„,e(J) X F4J) x F,,(J)"+i x Fe(J)"+2, 

where 

Fn,e(J) = ip(J X nr+\ F,(J) = iJpi(J;iJ-i(f^)) nLp(J;iJpi(f^)), 

and 

Fe(J) = T4-pi/2-i/2p( J; i^(s)) n Lp{J; W^-'/p{J:)), 

F„(J) = W^-'/'P{J; Lp(E)) n Lp( J; M^^^/p(E)). 

Then by locahzation and coordinate transformations it fohows from the maximal 
regularity result in [23] that the operator L defined by the left hand side of (|4.ip , 
(|4.2p , ()4.3p is an isomorphism from E into F x E^; see also Section 7. The range 
of L is determined by the natural compatibility conditions. If the time derivatives 
dt are replaced by 9t + w, w > sufficiently large, then this result is also true for 
J = M+. 

2. We introduce a functional analytic setting as follows. Set 

Xo - Lp,4n) X LpiQ) X w^'^/p{i:), 

where the subscript a means solcnoidal, and define the operator L by 

L(u, 7?, /i) = ( - {^i,/p)Au + V^p, -{d,/pK,)Ai!),-lpu ■ ly^j/lpj 
To define the domain D{L) of L, we set 

Xi = {(u,?9,/i)Gij2(f]\I])"+i x%3-i/P(s) :divu = in f^ \ S, 
-PsM = 1^1 = on E, M = a^?9 = on dQ}, 

and 

i?(L) = {(u,i?,/i) e Xi : PEb*^(u)j^E] = 0, {h/e,)j + ld,d,^^j = on E}. 
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Here j is given by j — [u ■ i^s]/|p^^l, and n is determined as the solution of the 
weak transmission problem 

(V7r|V0/p)2 = {{p,/p)Au\V(^)2, G H'p,in), 
lTrj = -<jA^h + 2lpi,{D{u)iy^\iy^)l on E, 
In/pj = 2l{pi,/p)iD{u)i^^\iy^)j - U^ on E. 

Let us introduce solution operators Tj, j e {1, 2, 3} as follows 

-Vtt = Ti{{fi,/p)Au) + T2{-aA^h + 2lp,{D{u)v^\v^)\) 
P 

+ T^{2l{pJp){D{u)v^\v^)}-Ud). 

We refer to Kohne, Priiss and Wilke [13] for the analysis of such transmission 
problems. 

Then the linearized problem can be rewritten as an abstract evolution problem 
in Xq. 

z + Lz = f, t>0, z{0) = zo, (4.4) 

where z = {u,'d,h), f = (fujejh), zo = (uQ,-&o,ho), provided ga = gu = 9 = 
99 = 9h — 0- The linearized problem has maximal Lp-regularity, hence (J4.4I) has 
this property as well. Therefore, by a result due to Hieber and Priiss, —L generates 
an analytic Co-semigroup in Xq; cf. Priiss [18], Proposition 1.1. 

Since the embedding Xi ■^ Xq is compact, the semigroup e~^* as well as the 
resolvent (A + L)^^ of —L are compact, too. Therefore the spectrum a{L) of L 
consists of countably many eigenvalues of finite algebraic multiplicity, and it is 
independent of p. 

3. We concentrate now on the case l^, ^ 0. Suppose that A with Re A > is an 
eigenvalue of — L. This means 

Xpu — p^Au + Vtt = in fJ \ E, 

divu = in f7 \ E, 

PsM = onE, 

^lp,D{u)iy^j + H + aA^hi^^ = on E, (4.5) 

u = on dH,, 

XpK.,-d-d.,Ad^Q inf}\E, 

|i?l =0 on E, 

^{h/e,)lu-j^^j/ll/pj~ld,d,^dj=Q onE, (4.6) 

d^i9 = on dn, 

l-p-^T{u,Tr,^)iy^-iy^l+h{l = onE, 

Xlpjh ~ Ipu ■ 1^^] = onE. (4.7) 
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Observe that on S we may write 

Uk= Pi2U + \hv^+ jv^/pk, fc = l,2. 

By this identity, taking the inner product of the problem for u with u and inte- 
grating by parts we get 

= A|pi/2y|2 _ ^^jy ^(-^^ ^^ ^^|^>)^ 

= A|p^/2y|2 ^ I T{u,TT,-d) : \/udx 

+ / (r2(-U2,7r2,^2)j^s •■"2 - ri(ui,7ri,i9i)z/s • Wi)dS 

- X\p"''u\l + 2|/iy'i?(,i)|2 + aA(^s/i|Ms + h{m^. 

since |T(u,7r, ??)i^5]] = aA^hv-^^ and |T(u, tt, i?)!^^ • ;/s/p| — h'd. On the other 
hand, the inner product of the equation for d with i? by an integration by parts 
and ld\ = leads to 

= A|(pK,)i/2^|2 -f \dl'^V-d\l + {{d^d^^dWd)^ 

where we employed {h/9^,)j = — [d* 9^^, ?9] . Adding theses identities and taking 
real parts yields the important relation 

= Re A|p^/^u|2 + 2\pl^^D{u)\l + aRe X{Ash\h)s 

+ e4Re\\{pK,y/^d\l + \dV^W^\l). (4.8) 

On the other hand, if /? :— ImA ^ 0, then taking imaginary parts separately we 
get with a ~ 1*{^\J)t, 

= I3\p^^^u\l - (Tl3iA^h\h)^ + Ima, 
= /3|(pK*)i/2^|2+Im a/61*, 

hence 

a{A^h\h)^^\p'/^u\j-9,\{pK,y/'d\l 
Inserting this identity into (|4.8p leads to 

= 2ReX\p^/^u\l + 2\py^D{u)\l + e^ldV^VdH 

This shows that if A is an eigenvalue of —L with ReA > then A is real. In 
fact, this identity implies ■(? = constant and D{u) = 0, hence j = by the Stefan 
condition, and then u = by Korn's inequality and the no-slip condition on dfl, 
as well as "i? = /i = by the equations for -d and h, since A 7^ 0. 
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4. Suppose now that A > is an eigenvalue of —L. Then we further have 
A / hdH = / (wfc • vy. - i/pk)dll = -pjT^ / jdS, 

hence the mean values h of h and j of j vanish since the densities are non-equal. 
Moreover, the identity 



implies also (z^lpK*)^^ = 0. Since Aj: is positive semidefinite on functions with 
mean zero in case S is connected, by (|4.8p we obtain u = -d = h = 0, i.e. in this 
case there are no positive eigenvalues. On the other hand, if S is disconnected, 
there is at least one positive eigenvalue. To prove this we need some preparations. 

5. First we consider the heat problem 

pK*A^-rf*Ai? = 0, xen\j:, 

-K9.,^]=5, xeE, (4.9) 

d^^ = 0, x€ an, 

to obtain -d — N^ g, where N^ denotes the Neumann-to-Dirichlet operator for this 
heat problem. The properties of N^ are summarized in the following proposition; 
see [55] for a proof. We denote by e the function which is identically one on E. 

Proposition 4.1. The Neumann-to-Dirichlet operator N^ for the diffusion prob- 
lem (|4.9p admits a compact self-adjoint extension to -/j2(S) which has the following 
properties. 
(i) If "d denotes the solution of ()4.9p . then 

in particular, N^ is injective for A > 0. 

(ii) For each a G (0, 1/2) and Aq > there is a constant C > such that 

iN^9\9)L,m > ^W^9\Uy), 9 e L^{^), A > Aq; 

hence 

C 
Wx le(L2(s)) - 1^' -^ - ^0- 

(iii) On i2,o(S) := {g G ^2(2) : ig\e)L^(^^^ = 0}, we even have 



iN^g\9)L,iy) > ^^;^|A^i'.9lL(E), 9 e i2,o(S), A > 0, 



c 

In particular, for A = 0, (j4.9p is solvable if and only if {g\e)L2i'S) — 0, and then 
the solution is unique up to a constant. 
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6. We also need a corresponding result for the Neumann-to-Dirichlet operator N^ 
for the Stokes problem. It is defined as follows. Given a function g e Wp'^ (S); 



(4.10) 



and define N^g :— u ■ vy: on S. For this well-defined operator we have 

Proposition 4.2. The Neumann-to-Dirichlet operator N^ for the Stokes problem 

(14.101) admits a compact self-adjoint extension to -/j2(S) which has the following 

properties. 

(i) If u denotes the solution of (j4.10p . then 

iM9\9)L,i^) = A / p\u\^dx + 2 f ti^,\D{u)\ldx, A > 0, .9 e i/^'lS). 
Jn Jn 

(ii) For each a G (0, 1/2) there is a constant C > such that 



solve the Stokes problem 




pXu — fi^Au + Vtt = 0, 


xefi\E, 


divu = 0, 


xen\^, 


N = 0, 


xeY., 


-lT{u,n,'d)uY,] ^ gvT., 


a; e S, 


u^Q, 


xedn, 



(iVf5lff)L.(E) > ^^y^" l^fglL(E). 9 e i2(S), A > 0. 



In particular, 



l^f kL.(s)) < (^-^, A>0. 



(iii) Let E^ denote the components of S and let e^ be the function which is one 
on Sfe, zero elsewhere. Then {N^g\ek)L2{S) = for each g e L2(S). In particular, 
N^Ck = for each k, and N^g has mean value zero for each g £ L2(S). 

This result is proved in [27], Proposition 4.1 for the case of equal densities. It 
carries over directly to the case fp] 7^ considered here. Note that N^ even on 
L2,o(S) is not injective in case S is disconnected. 

7. Next we solve the asymmetric Stokes problem 

pXu — /i* Au + Vtt = 0, a; e fi \ S, 

divu = 0, x £ fl\'E, 

P^M^P^lTiu.n^i^y^j^O, xeE, (4.11) 

-|r(u,7r,'i9)i/s • J'e] =5i, a; e s, 

-|T(u,7r,^)i^s • i^s/p] =52, a; € E, 
to obtain as output 

k^lpu- z^sl/M = ^^51 + ^1^52, J^lu- ^e]/[1/p1 = ^f 51 + ^f 52. 
For this problem we have 
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Proposition 4.3. The operator Sx for the Stokes problem ()4.11|) admits a bounded 
extension to L2.q{T,)'^ for A > and has the following properties. 
(i) If u denotes the solution of (j4.1ip . then 

{Sx9\9)l,(Sp = a / p\u\^dx + 2 f fi,\D{u)\ldx, A > 0, .g G L^A^f r)Hl^^{j:)\ 
Jn Jn 



(ii) Sx G B{L2fi{^) ) is self-adjoint, positive semidefinite, and compact; in par- 
ticular 

sv-[sVY, sf^[sfr, s'^^[sfr. 

(iii) S^ and S^ are injective in L2,o(S), and with Gx — [S'^'^]^^ we have 

N^ = S^ - S^-G\S^ . 

Gx is self-adjoint and positive definite on L2,o(S), its resolvent is compact in 

L2,o(S); for each A > 0. 

(iv) For each /? G (0, 1/2) there is a constant Cp > such that 

\Sx\b(L2,o(^Y) < {ij^x)P ' ^-'^' 

(v) l'5'A|e(L2,o(s)2,_f/i(s)=^) ^ ^ uniformly for A > 0. 

(vi) S]^ , Sj^ : 2^2,0 (S) — >■ -^2 (^) ^ i2.o(S) are isomorphisms, for each A > 0. 

Proof, (a) First observe that for the traces Uj of u on S we have 

Uj = Pem + fci^s + JVT,! P] = Mb + JVT.1 P] on S. 

To prove assertion (i), let (u,7r) denote the solution of (|4.1ip . Multiply with u 
and integrate by parts to the result 

A / p\u\ dx -\- 2 ii^,\D{u)\2dx — / div {T{u,tt, {>)€,) dx 
Jn Jn Jn 

= - flu- T{u, t:, ^)Md^ = - f IK+jVs/p) • T{u, n, ^)i^j:]dJ: 
uE ■ lT{u, n, z9)i/sldS - / ]lT{u, n, §)v^ ■ v^/ p}di: 

= / kgidT,+ / jg2d^^ {g\Sxg)L2iT,y- 

A similar computation yields 

('5'a5|^)l2(s)2 = (.9|'S'a/i)l2(s)2, 

hence Sx is self-adjoint in i2(S)^, thereby proving the first part of (ii). 

(b) To prove injectivity of 5'|^, let gi = and j — 0. Then (i) implies Xu = 
D{u) — hence Vtt = in fi \ E, and so tt is constant in the components of the 
phases. Next 

= ,91 - - lT{u, TT, ^)u^ ■ ly^j = H , .92 - - lT{u, TT, ^)u^ ■ u^/pj - Itt/pJ 
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shows that n is even constant in the phases, and so (72 is constant on S hence zero 
since its mean value vanishes. Injectivity of Sj^^ is shown in a similar way. This 
proves the first assertion in (iii), the second one follows from the definitions of N^ 
and S\, and the third one is a consequence of (ii). 

(c) To establish the boundedness properties of Sx, we proceed as follows. Suppose 

1/2 
51, 32 e £2,0 (S) n H2 (S) are given. We decompose the pressure into n = q + qo, 

where the one-sided traces Qq are uniquely determined on E by 

M=gi, lqa/pl=92; 

then q^ e L2fi{T,). Extend go to all of O in iJj (^ \ 5]), by solving the problem 

Aw = in O \ E, v — qo on Y,, w = on dil. 
There is a constant co > such that 

l'?oli/i/2(fi\s) ^ co|5|l2(s)- 

Then (u, q) satisfies the asymmetric Stokes problem 

pXu — fi^,Au + Vq = —Vqo, a; e 51 \ S, 

div u — 0, X G fi \ S, 

Pj:luj=Pj:lfi,D{u)i^^j^O, a; e S, (4.12) 

~2lfi,D{u)uj: ■ i^s] + M = 0, xe S, 

-2l^i,D{u)iy^-jy^/p] + lq/p]=0, a;eS, 

M = 0, X E dft. 

Note that the interface conditions are now homogeneous. Let —A denote the gen- 
erator of the associated analytic Co-semigroup in L2,cr(51), which is exponentially 
stable. Then we have 

ux:=-{\ + A)-'[Vqo/p], 

and with 

Sxg = (fc,j) = ilpux ■ /.sl/M, lux ■ ^sl/Il/p]) 

we estimate using trace theory and the resolvent estimate for A and Re (t{—A) < 

as well as 

I^a5|l.(e) < C|«AUj(fAE) < C{1 + A)^-i|go|^i/2(,,^5,) < C{1 + A)^-i|.gU,(s), 

with (3 — {a + l/2)/2, for each fixed a > 1/2. These estimates are uniform w.r.t. 
A> 0. 



24 J. PRUSS, S. SHIMIZU, AND M. WILKE 

(d) The proof of (vi) is more involved. To obtain surjectivity of S\^ we have to 
solve the problem 

pXu — /i* Au + Vtt = 0, a; e ri \ S, 

divu = 0, X £ fl\T,, 

P^lul = P^lT{u,n,i9)iy,:l = 0, x G S, (4.13) 

Ipu ■ j/sl = [pp, a; e E, 

-|T(u,7r, i?)z^E • i^s/pl = 0, a; e S, 

u — 0, X £ dfl, 

for given k e -ffl(S) n L2^o(S) with output [Sl^r'^k = gi = -|r(M, 7r,^)i/E • J^sl- 
Similarly, to prove surjectivity of S^^ the problem to solve is 

pXu — /i* Am + Vtt = 0, a; G fi \ S, 

divu = 0, a; G ri \ S, 

P^lul = P^lT{u,n,i})iy^l = 0, a; G S, (4.14) 

|w • i^s] = [l/p]j a;GS, 

-|r(w,7r,^)i.s-J^sl-0, a;GE, 

u = 0, a; G 9S1, 

for given j and the output will be [^I^J^^j ~ g2 ~ ~[r(u,7r, z?)j/s • i't./pI- These 
problems can be solved in a similar way as in the proof of Theorem 13.11 Using 
the method of localization and perturbation we reduce the problems to the case 
of a flat interface, i.e. ft — R", E — R"^^ x {0}. Then we may employ a partial 
Fourier-transform in the same way as in the proof of |23) . Theorem 3.1. D 

8. Now suppose that A > is an eigenvalue of L. We set g = l^,j/9^, , gi — —aA^h 
and 52 — ~l*'d = —c^N^j, c* — l1/6^ > to obtain 

3 = -SfaA^h - c^SfN^j, Xh - -Sl^aA^h - S^^c^N^j. 

Observing that / + Ci,S^N^ is injective, hence boundedly invertible in L2,o(S) 
by compactness of S\ , we may solve the first equation for j to the result 

3^^a{I + c,SfN^)-^SfA^h. 

Inserting this into the second, the equation for h becomes 

= A/i + (7(511 - Sfc,N^[I + c,SfN^)-^Sf)A^h, 

or equivalently with R\ = G\S^ 

Q = Xh + a{N^ + Rl{c,N^ + Gx)-^Rx)A^h. 

Next we observe that the operators N^ + R*^{c^,N^ + G\)^^R\ are injective for 
A > 0; in fact if 

N^h + Rlic^N^ + Gx)-^Rxh = 0, 



INCOMPRESSIBLE TWO-PHASE FLOWS WITH PHASE TRANSITIONS 25 

then with v\ = {c*N^ + G\)^^R\h, formmg the mner product with h in L2,o(S) 
we obtain 

= {N^h\h)j: + {Rlic^N^ + Gx)-'Rxh\h)^ 

> c\Nih\l + {{c,N^ + Gx)vx\vxh 

hence N^h = and v\ — 0, which imphes Sf^h = S\^h = 0, hence h — hy 
injectivity of S^^. Setting now 

Tx = [N^ + R*xic*N^ + Gxr'Rxr\ 
we arrive at the equation 

XTxh + aA^h^O. (4.15) 

A > is an eigenvalue of — L if and only if the problem (|4.15p admits a nontrivial 
solution, i.e. if and only if is an eigenvalue for Bx :— XTx + oAy.. Here the 
domain of Bx is that of Ay, , Ta is a relatively compact perturbation of Ayi ■ 
We consider this problem in ^2.0 (5^)- Then As is selfadjoint and 

<j{Ash\h)s > -^^^\h\l. 

On the other hand, N^ , N^ are self-adjoint, compact and positive semidefinite on 
L2,o(S), and Gx is self-adjoint, positive definite, and has compact inverse. Hence 
Tx is self-adjoint, positive semidefinite and T^ is compact as well. If /i > is an 
eigenvalue of Tx then 

hence we get 

^i-'\hk<c\h\s.. 

since by Propositions 14.31 and 14.11 ^a, N^ and (/ + c^,Sl^N^)^^ are bounded in 
L2,o(S), uniformly for large A. Therefore fi — //(A) > co > 0, for large A, and so 

(j(n~l) 



{Bxh\h)s = X{Txh\h)s + cr{Ash\h)s > (cqA - ^ ^, ' )\h\ 



IS- 



This proves that Bx is positive definite, hence (J4.15I) has no nontrivial solution, 
for large A. 

But for small A > we have with h = X]fc=i XSt^fe, where XA denotes charac- 
teristic function of a set A, hk = const, on S^, X]fe=i ^k — 0, 

KTxh\h) ~ ^(!l^^„i?"-i 5] /i^ < 0, 

k 

since with h £ H^^^{T.) n L2,o(S) 

Txh={l-c4Sir'Sl'N^{I + c.SrN^r'Sir' 

■.= {I-Kx)-'[SiY'h^Toh 

in 2^2,0 (S) as A -T> 0. This can seen from the mapping and continuity properties of 
Sx and N^; Kx is compact and continuous in S(L2,o(S)), and / — Kiy is injective, 
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hence invertible and its inverse is continuous as well. This shows that Bx has a 
nontrivial kernel for some Ao > 0, which implies that — L has a positive eigenvalue. 
Even more is true. We have seen that B\ is positive definite for large A and 
Bq = cfAy. has —a{n — 1)/ R^ as an eigenvalue of multiplicity ?7i — 1 in L2,o(S). 
Therefore, as A increases to infinity, m — 1 eigenvalues fikW of B\ must cross 
through zero, this way inducing m — 1 positive eigenvalues of ~L. 

9. Next we look at the eigenvalue A = 0. Then (|4.8p yields 

hence i) is constant, D{u) = and j = by the flux condition for i). This further 
implies that |u] = 0, and therefore Korn's inequality yields Vm — and then we 
have u = by the no-slip condition on dQ. This implies further that the pressures 
are constant in the phases and [tt] = —crAsh, as well as l^^ = — |7r/p]. Thus the 
dimension of the eigenspace for eigenvalue A = is the same as the dimension of 
the manifold of equilibria, namely m-n+2 if ili has m > 1 components. The kernel 
of L is spanned by eg — (0, 1, 0), e/i = (0, 0, 1), ejk — (0, 0, Yj') with the spherical 
harmonics Y^^ of degree one for the spheres T,'^ , j = 1, . . . ,n, k = 1, . . . ,m. 

To show that the equilibria are normally stable or hyperbolic, it remains to 
prove that A = is semi-simple. So suppose we have a solution of L{u,'d,h) = 
^j.k c^jkSjk + (3e0 + jch. This means 

div u — 0, X £ fl\T,, 

M-lp-ilji^s, xeJ:, (4.16) 

M = 0, X & dfl. 

-d*A^ = /3, xen\^, 

1^91 =0, X e E, 
-{h/e,)j - |49,^z91 =0, X e S, (4.17) 

d^i) = 0, x e on. 

l-p-^T{u,TT,d)iy^-iy^j + hd^O, a; G S, (4.18) 

-Ipu ■ ly^j/lp] - J2 "^■^' + ^' i > 0, a; e E. 

j,k 

We have to show ajk = /3 = 7 = for all j. Integrating the divergence equation 
for u over H. we get J^jdT, = 0, and integrating that for -d yields 

P\n\ = - /" d^AMx ^ f Id^d^^^dT. ^ -{h/e,) f jdT. = 0, 
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which implies /3 = 0. As in 3. we then obtain 

2\,il'^D{u)\l + {lp-^T{u,n,d)v^ ■ M\j)^ = 0, 
and 

Adding these equations yields 

This implies D{u) — 0, d constant, hence j — since l^ ^ 0, and then u = by 
Korn's inequality, which in turn yields 

Thus 7 — ajk — for all j, k since the spherical harmonics Yj^ are independent 
and have mean zero. Therefore eigenvalue A = is semi-simple. 

10. Let us consider now the case l^ — 0. Then the temperature equation 
decouples completely from that for u and h. It only induces one dimension in the 
kernel of L, but no positive eigenvalues. In this case, as now c* = the derivation 
in 8. yields the equivalent problem 

Xh + (7Sl^AY:h = 0. 

As S\^ is positive semidefinite and injective, this equation admits no nontrivial 
solutions if Re A > and A is non-real. If A > then Tx = [S\^]~^ and may employ 
the same arguments as in 8. to obtain the same conclusions as in case I* 7^ 0. 
Finally, for A = we have the same kernel of L as before, and is semi-simple for 
L, too. 

11. Let us summarize what we have proved. 

Theorem 4.4. Let L denote the linearization at e* :— (0,0*, E) ^ £ as defined 
above. Then —L generates a compact analytic Ca-semigroup in Xq which has 
maximal Lp -regularity. The spectrum of L consists only of eigenvalues of finite 
algebraic multiplicity. Moreover, the following assertions are valid. 

(i) The operator —L has no eigenvalues X ^ with nonnegative real part if and 

only if E is connected. 

(ii) // E is disconnected, then —L has precisely m — 1 positive eigenvalues. 

(iii) X = is an eigenvalue of L and it is semi-simple. 

(iv) The kernel N{L) of L is isomorphic to the tangent space Te,£ of the manifold 

of equilibria £ at e^. 

Consequently, e* — (0,0*, E) G £ is normally stable if and only if T, is connected, 
and normally hyperbolic if and only if E is disconnected. 
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5. Nonlinear Stability of Equilibria 

1. We look at Problem (jl.l|) . (jl.2p . (jl.3|) in the neighborhood of a non-degenerate 
equilibrium e^ — (0,0*, F.^) S i5. Performing a Hanzawa transform with reference 
manifold E = P.^ as in Section 3, the transformed problem becomes 

pdtu — /i*A?i + Vtt ~ pFu{u, ??, ft, tt), in J7 \ S, 

drv u — Gd{u,h), in J7 \ S, 

PsM=G„(u,z9,/i), onE, 

-PE[M*(Vw+[Vu]T)li.s-G,(w,^,/i), onE, (5.1) 

-{T{u, TT, i9)i^s • M + cr^s/i = Gu{u, i9, ft) + G^(ft), on E, 

/,^ - lr(w, ^, i9)z.E • i^s/p] = Gh{u, z9, ft), on E, 

M = 0, on 9f2, 

w(0) = Wo, in ri, 

where /^* = /x(6'*) > and At. = -H'{0) = -{n - l)/Rl - As. 
For the relative temperature d ~ {9 — 9^,)/6^, we obtain 

pn^dt^-d^M^ Fg{u,d,h), in O \ E, 

m = 0, on E, 

-{k/9,)lu-i,TJ/ll/pj-ld.d,^n^Ge{n,^,h), on E, (5.2) 

a^i? = 0, on dn, 

i9(0) = i?o, in n, 

with K* = k{9^,) > and Z* = /(fi**) 7^ 0. Finally, the evolution of ft is determined 

by 

dth-lpu-iyTl/lpl^Fh{u,h), onE, (5.3) 

ft(0) = fto. 

Here the nonlinearities are of class C^ from E to F, and satisfy ^'(0) — G^(0) — 
for all j — u,9,h and k — d,u,T^iy,^^h,9. Let w :— {z,t:) :— {u,'d,h,T:) and 
zo :— (uo, 1^0, ^o). We will frequently make use of the shorter notation 

Lw = N{w), z(0) = Zo 

for the system ((^ - (|0)) . 

The state manifold locally near the equilibrium e* — (0, 0*,P*) reads as 

SM :={(M,79,ft) e Lp(r!)"+i x ^^(E) : {u,d) e W^-^^P{n\j:y'+\ h e W^-^/p, 

^sM =G„(u,^,ft), -Fs[M*(VM+[Vu]T)lzys = G,(M,i?,ft) on E, 
divu = Gd(u,ft) in 17 \E, u = 9^61 = on 917, (5.4) 

1^1 = 0, -{U/9,)lu-,.T]/ll/pj - Id*a,,^l = Geiu,d,h) on e}. 
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Note that due to the compatibihty conditions this is a nonhnear manifold. We 
need to parameterize this manifold over its tangent space 

SX := [{u,d,h) e £p(0)"+i X C2(E) : (u,^) E W^-^^P{n\j:)"+\ HeW^'^'p, 

PsM - 0, -Ps[m*(Vm + [VuY)\v^ = on S, 

div u = in f7 \ S, u = 9^61 = on dVl, (5.5) 

m = 0, -{l*/e*)lu ■ M/1^/P\ - ld*d..n = on s}. 
We mention that the norm in SX is given by 

\{u,d,h)\sX = |u|^2-2/p + |^|^2-2/p + |/l|^3-2/p. 

2. To parameterize the state manifold SAi over SX near the given equilibrimn 
we consider the linear elliptic problem 

puju ~ fi^Au + Vtt = in f2 \ S, 

div u ~ g^ in f2 \ S, 

pK,wi9-d,A^ = mn\T,, 

u = d^d = on dfl, 

PsM=5« onS, (5.6) 

-PslAi*(Vu + [Vu]'^)P^ = gr on E, 

|?9] =0 on E, 

-{h/9.)lu ■ ^e]/[1/p1 - Id.d^^n = 5fl on E, 

-|r(M,7r,'i9)i/E • i^s] = 5^ on S, 

/*i9-|T(u,^,z9);/s-i^s/p]=0 on S, 

for given data gd, gu, dr, 9^, ge- For this problem we have the following result. 

Proposition 5.1. Suppose p > 3 and letuj > be sufficiently large. Then problem 
(|5.6p admits a unique solution (u,7r,^,ft,) with regularity 



if and only if the data {gd, gu, gr, gu, ge) satisfy gu ■ v^ = gr ■ i^^ = and 

gd e w;,-'^P{n\^)nH;\n), {gr,g.,ge) e %i-=^/^(E)"+^ g„ e ^^-^^/^(S]). 

T/ie solution map {fd, gu, gr, 9^, ge) *— >■ (i^iTTj^?) is continuous in the corresponding 
spaces. 

This purely elliptic problem can be solved in the same way as the corresponding 
linear parabolic problems. 
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3. For the parametrization, fix any large a; > 0. Given (u, ■&, h) e SX sufficiently 
small, and setting (u, •&, h) = (m, i9, h) + (u, d, 0), solve the auxiliary problem 

puiu — /i* Au + Vtt = 0, in rj \ E, 

drv u — Gd{u,h), in i7 \ E, 

pK^,ujd — di.A^ = 0, in i7 \ E, 

u = d„d = 0, on 951, 

-PeN=G„(u,^,/i), onE, 

^PslAi*(VS+[Vfi]T)]z.s-G,(w,^,/i), onE, 

m =0 on E 

~{L/e,)lu ■ ly^j/U/pj - ld,d,Jj = Ge{u,^,h), on E, 

-|T(M,7f, i?)z^E • i/sl = -o-^e/i, on E, 

hd~lT{u,Tt,^)iy^-i^s/pj^O, onE, 

by means of the implicit function theorem, employing Proposition 15. II Then with 
z = (u, ??, 0) and z = {u,'d,h), z = z + z, we obtain z = 0(5) , with a C-'^-function 
such that 0(0) = 0'(O) = 0. Then z ^ z + 0(z) e 5A^, hence SM is locally near 
parameterized over SX . To prove surjectivity of this map, given (u, 1}, h) G SAi 
sufficiently small, solve the linear problem 

puju — p^Au + Vtt = 0, in 57 \ E, 

div u — Gd{u,h), in 57 \ E, 

pKi,uj-d — d^A-d = 0, in 57 \ E, 

u = d^-d — 0, on 957, 

-PeN=G„(u,^,/i), onE, 

^P^{p^,{Vu+[VuY)\u^ ^ Gr{u,d,h), onE, 

m =0 on E 

~{h/B,)lu ■ M/{l/p} ~ {d,d,J\ = Ge{u,^,h), on E, 

-|T(w,7f, ??)z^s ■ ;^sl = -cr^s/i, onE, 

hd-lT{u,Tr,d),yj:-i^j:/pj=0, onE, 

and set z — z — z. Then we see that z = (^{z), hence the map z i— > z + (j){z) is also 
surjective near 0. 

4. Next, let Wao ■= (zoc^'n'oc) ■= (^oo, 'Joo, ^ocTToo) be an equilibrium of (IS.ip - 
(j5.3p such that Zoo — z^o + 4>{zco) G SM is close to the part z* = (0, 0, 0) of the 
fixed equilibrium w^ — (0, 0, 0, tt*). Clearly, Uao = 0, t?oo = const., tt*, tToo = const, 
and it holds that Luiqo = N{wrxi)- 
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Given a solution w e E of Problem (|5.ip - (|5.3p . we decompose w as w = w + 
w + Woo , where now w solves the following auxiliary problem 

ujpu + pdtu — n^Au + Vtt — pFu{u,i}, h, tt), in 51 \ E, 

div u = Gd{u, h), in fl \ S, 

pn^uj^ + pn^dt^ - d^A^ = Fe{u,-d,h), in rj \ S, 

u — d^d — 0, on (90, 

Py.lu}^Gu{uJ,h), onE, 

(5.7) 

-P^[p,{Vu+[VuY)}v^ = Gr{u,d,h), on S, 

|t91 = on S 

~{h/e,){u ■ MII^IpI - {d*d.J\ = Ge{u, d, h), on S, 

-lT{u, Tt, i9)i^s • i^sl + <jAj:h = G^{u, i9, h) + G^{h) ~ G^{h,^), on E, 

hd - lT{u, 7f, z9)iys • i^^/p] = Ghiu, i9, h), on E, 

w/i + St/i- |pu- i^sJ/M =i^/i(w, /i), on E, 

u{o) = uo, m^^o, m^ho, I 

where zq = zq + zq = zq + (f>{zQ). The remaining problem for w reads 

pdtu — /i* Am + Vtt = ujpu, in fl \ S, 

divu = 0, in fl \ E, 

pn^.dt'd — d^AiS = pK^.ui'd, in fl \ E, 

u = di,^ = 0, on 9ri, 

PeM = 0, -PT.l^i.{Vu + [VuY)\y^ = 0, on E, (5.8) 

1^1 =0 on E 

{hie.)lu ■ M/I^/Pl - {d*d.J\ =0, on E 

-lT{u,^,d)vT.-M + <^A^h = 0, onE, 

lJ-lT{u,^,d)v^-v^/p} = 0, onE, 

dth-\pu- v^}/\p\^ bjJi, onE, 

u(0) = Mo, ^(0) = ?9o, /J(0) = /lo- 

5. Problem (|5.7p can be written abstractly as 

l^^w = N{w + w + Woo)- N{woo)- (5.9) 

And (|5.8p abstractly becomes the evolution equation 

S + LS^R{z), t>0, z(0) = zo, (5.10) 
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in the Banach space Xq, with L defined as in Section 4 and R{z) :~ lo{{I — 
Ti)u, K^,'3,h). This abstract problem can be treated in the same way as in the 
proof of Theorem 5.2 in Priiss, Simonett and Zacher [2S]. This impUes the following 
result. 

Theorem 5.2. Let p > n + 2, a,pi,p2 > 0, pi ^ p2, and suppose ipj G ^^"^(0,00), 
Pj,dj G C^(0,oo) are such that 

Kj{s) ^ -sipj{s) > 0, pj{s)>0, dj{s)>0, s € {0,00), j ^1,2. 

Then in the topology of the state manifold SA4 we have: 

(i) (0, 6'*,r*) G E is stable if and only ifT^, is connected. 

(a) Any solution starting in a neighborhood of a stable equilibrium exists globally 

and converges to a probably different stable equilibrium in the topology of SA4. 

(Hi) Any solution starting and staying in a neighborhood of an unstable equilibrium 

exists globally and converges to a probably different unstable equilibrium in the 

topology ofSAd. 

6. Global Existence and Convergence 
In this section we study the global properties of problem (|l.ip . (|1.2p . (|1.3p 

6.1. The Local Semiflow. We follow here the approach introduced in Kohne, 
Priiss and Wilke [T3] for the isothermal incompressible two-phase Navier-Stokes 
problem without phase transitions and in Priiss, Simonett and Zacher [26j for the 
Stefan problem with surface tension. 

Recall that the closed C^-hyper-surfaces contained in fl form a C^-manifold, 
which we denote by A4'H^{^). The charts are the parameterizations over a given 
hyper-surface E according to Section 2, and the tangent space consists of the 
normal vector fields on E. We define a metric on A^'H^(ri) by means of 

d^«2(Ei,E2) ■.= dH{N^^l,N^^2), 

where dn denotes the Haussdorff metric on the compact subsets of M" introduced 
in Section 2. This way Ai'H^{^) becomes a Banach manifold of class C^. 

Let d^(x) denote the signed distance for E.We may then define the level function 
ifs by means of 

(^s(a;)=0(ds(x)), x e R", 
where 

(j){s) — s(l — x('5/a)) + sgnsxis/a), s G K. 

It is easy to see that E = ip^^{0), and \7(pY:ix) — i^T:ix), for each a; e E. Moreover, 
K = is an eigenvalue of V^(/3s(a;), the remaining eigenvalues of 'S/'^cpsix) are the 
principal curvatures Kj of E at a; € E. 

If we consider the subset M'H^{fl,r) of A^H^(r2) which consists of all closed 
hyper-surfaces F G A^'H^(r2) such that F C f2 satisfies the ball condition with 
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fixed radius r > then the map $ : Mn'^{n,r) -^ C'^iCl) defined by $(r) = ipr 
is an isomorphism of the metric space ^^7^^(17, r) onto $(A^'H^(il,r)) C C^{ri). 

Lets-{n-l)/p>2;ioTr eMn^{n,r),wede&neT e W^{n,r)iiipr G W^{n). 
In this case the local charts for F can be chosen of class W^ as well. A subset 
A C Wp{fl,r) is (relatively) compact, if and only if $(A) C Wp{fl) is (relatively) 
compact. 

As an ambient space for the state manifold SAi of Problem (|l.ip . (|1.2p . (I1.3P 
we consider the product space Lp(J7)"+^ x A^H . 

We define the state manifold for the problem SAi as follows. 

SM := [iu,e,T) e Lp{Qy'+^ x Mn^i^) : {u,e) e w^-^^pin\Ty'+\ r e Typ3-2/p^ 

divu^O in f7, 9 > in Q, u ^ d^9 ^ on 9fi, (6.1) 

10] = |Prw] = Pr[/i(e)(VM + [Vuf)liyr =0 on T, 
/(^)|u • H/[l/p] + K0)a.,01 = on r}, 

Charts for these manifolds are obtained by the charts induced by AiJi (51), fol- 
lowed by a Hanzawa transformation. 

Applying Theorem 13.21 and re-parameterizing the interface repeatedly, we see 
that (fLT|) . ([L2]), ([0|) yields a local semiflow on SM. 

Theorem 6.1. Let p > n + 2, a,pi,p2 > 0, pi "^ Pi, and suppose ipj G ^^"^(0,00), 
Pj,dj G C^(0,oo) such that 

Kj{s) = -sipj{s) > 0, pj{s)>0, dj{s)>0, s G (0,00), j = 1,2. 

Then problem lll.l]) . \1.2\) . Ill.3\) generates a local semiflow on the state manifold 
SM. Each solution {u,9,T) of the problem exists on a maximal time interval 
[0,t*), whereto = t*(uo, ^0, To). 

Note that the pressure tt as well as the phase flux j are dummy variables which 
are determined for each t by the principal variables (u, 9, T). In fact, j is given by 

j^lu-i^rVll/pl 
and TT is determined by the weak transmission problem 

(V7r|V0/p)i.,(o) = i2p-'divipi9)Diu)) - u ■ Vti|V0)i,(a), G Hl,{n), 

M = 2lp{9)D{uyr ■ M + aHr - U/pjf, 

In/pj = 2l{p{9)/p)D{u)iyr ■ z^rl - ll/2p'lf - [V^^)! on F, 

where D{u) — (Vu + [Vu] )/2 as before. Concerning such transmission problems 
we refer to |Ti]. 
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6.2. Compactness of Orbits. Let us briefly comment on tlie fact of liow one 
can use tiie time weights t^~^, jjl £ (1/p, 1] introduced in Subsection 13.31 to prove 
relative compactness of bounded orbits in SA4. For simplicity we assume that 
the free boundary is fixed, hence we do only consider the orbits {u(t)}tgro^t,) and 
{^(^)}te[o.t*) of tti6 velocity field u and the temperature field 6, respectively. 
Assume that there exists M > such that \u{t)\^2~2/p + \6{t)\^2-2/p < M for 

all t G [0, t*). Since n is bounded it follows that {M(0}tG[o,t.) and {^(i)}te[o,t.) are 
relatively compact w.r.t. the topology of Wp'^~ . The continuous dependence of 
the solution on the initial data in Wp'^ and the instantaneous regularization of 
the solution, see Corollary 13. 3| yield that the orbits of u and 9 are also relatively 
compact in Wp . This in turn yields that the solution exists globally, i.e. 
t^, — oo. We refer to [13] for more details. 

If one considers free boundary problems, then one has to work in the setting 
which has been introduced in Subsection 16. II In particular, to show relative com- 
pactness of the orbits in SA4 one has to perform several Hanzawa transforms for 
each of the finitely many balls which cover the relatively compact set {r(t)}t>o in 
Wp~ ~'^{fl,r). The pull backs of the velocity and temperature field for each of 
these balls will then be relatively compact in Wp ". More details are given in 



the proof of Theorem 16.21 below. 

6.3. Convergence. There are several obstructions against global existence: 

• regularity: the norms of either u(t), 9{t), or T(t) become unbounded; 

• geometry: the topology of the interface changes; 
or the interface touches the boundary of J7. 

• well-posedness: the temperature becomes 0. 

Note that the compatibility conditions, 

div u{t) = in r? \ r(t), u = 9^6* = on dn, 

Pr[u(i)l = m = P^l^iDiu)it)] - /(0)|u • ^rVll/p] + ^6)8^,6] - on r(i), 

are preserved by the semifiow. 

Let (u, 9, F) be a solution in the state manifold SA4 with maximal interval 
[0,t*). By the uniform ball condition we mean the existence of a radius ro > 
such that for each t, at each point x G F(i) there exists centers Xi G ^i{t) such 
that Brg{xi) C fli and F(t) n Bro{xi) = {x}, i ~ 1,2. Note that this condition 
bounds the curvature of F(<), prevents parts of it to touch the outer boundary 
dft, and to undergo topological changes. Hence if this condition holds, then the 
volumes of the phases are preserved. 

With this property, combining the local semifiow for ()l.ip . ()1.2p . ()1.3p with the 
Ljapunov functional and compactness we obtain the following result. 
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Theorem 6.2. Let p > n + 2, a, pi, p2 > 0, pi ^ p2, and suppose ijjj G (^"^(O, oo), 
Pj,dj e C^(0, oo) are such that 

Kj{s) ^ -sip'-{s) > 0, pj{s)>0, dj{s)>0, s e {0,oo), j = 1,2. 

Suppose that {u, 6, T) is a solution of il.l\) . il.^) . il.S\) in the state manifold SM. 

on its maximal time interval [0,t*). Assume there is constant M > such that 

the following conditions hold on [0,i*).' 

(V W{t)\yy2-2/p, \9{t)\^2-2/p, \r{t)\^3-2/p < M < oo; 

(ii) r(t) satisfies the uniform ball condition. 

(iii) 0{t) >1/M onCt. 

Then i* = oo, i.e. the solution exists globally, and its limit set uij^.(u,6,T) C £ 
is non-empty. If further {0,Ooa,Too) G '-^+{'a,0,T) with Too connected, then the 
solution converges in SA4 to this equilibrium. 

Conversely, if {u{t),0{t),T(t)) is a global solution in SM which converges to 
an equilibrium {0,9.^,,T^) € £ in SM as t ^ oo, then (i), (ii) and (iii) are valid. 

Proof Assume that (i), (ii) and (iii) are valid. Then r([0,i*)) C Wp~'^'^{^,r) is 
bounded, hence relatively compact in Wp '^{Vl,r). Thus we may cover this 

set by finitely many balls with centers Efc real analytic in such a way that 

dist^y3-2/p-e(r(i),i;j) < 5, 

for some j = j{t), t G [0,t*). Let Jj. = {i e [0,t*) : j{t) = k}; using for each k a 
Hanzawa-transformation S^, we see that the pull backs {{u{t, ■), 6{t, •)) o S^ : i G 
Jfe} are bounded in Wp~ '^(r2\Efc)"+^, hence relatively compact in Wp~ '^^^{Vt\ 
'^k)^'^^ ■ Employing now Corollary 13.31 we obtain solutions {u^,9^ ,T^) with initial 
configurations {u{t),9{t),T{t)) in the state manifold on a common time interval 
say (0,t], and by uniqueness we have 

{u\T),e\T),T\T)) = {u{t + T),e{t + T),T{t + T)). 

Continuous dependence implies that the orbit of the solution {u{-),9{-),T{-)) is 
relative compact in SM, in particular i* = oo and [u, 9, r)(M_|_) C SM is relatively 
compact. The negative total entropy is a strict Ljapunov functional, hence the 
limit set lu+{u, 9, T) C SM of a solution is contained in the set £ of equilibria. By 
compactness lo+{u,9,T) C SM is non-empty, hence the solution comes close to 
£, and stays there. Then we may apply the convergence result Theorem l5.2l The 
converse statement follows by a compactness argument. D 

7. Proof of Theorems 3.1 

In this section, we prove Theorem l3.1[ employing the methods introduced in [2] . 
Actually, most of the arguments given in that paper remain valid for the problem 
under consideration here, hence we restrict on the necessary modifications. 



36 



J. PRUSS, S. SHIMIZU, AND M. WILKE 



7.1. Flat interface. In this subsection, we consider first the hnear problem with 
constant coefficients for a flat interface. Due to the jump in the velocity, this 
problem differs from that in ('Hj . 



dth 



pdtu - fioAu + Vtt = pfu 

div u = cjd 

{v\ + coVxh = Ps5« 

-lp,adyvj - l^oVj^wl = gv 

-2[iJ.odywl + [tt] - aAxh = g^ 

-2{pQdywlpl + {-k/pI = gh 

^bo-V^h/lp\^h 
m(0) = Mo 
h{Q) = /lo 



m 


K", 




in 


M", 




on 


R"- 


-1 


on 


M"- 


-1 


on 


K"- 


-1 


on 


R"- 


-1 


on 


M"- 


-1 


in 


R", 




on 


R"- 


-1 



(7.1) 



with ^0 > 0, Co e R, 6o G 



Here we have identified 



1 ^ M"-i X {0}, 

and R" — R" \ R"^-'^. It is convenient to split u — {v,w), g — {gv,gw) into 
tangential and normal components. Since the 3rd and the 7th equation differ 
from the coupled system (3.1) and (3.3) in [23], we solve (j7.ip with all right hand 
sides except fh- Employing Laplace transform in t and Fourier transform in the 
tangential variable x E R"^^ similar to Sec. 4 of [13], the transformed equation for 
h reads 



Xh- 



— fh- 



lpw{0)j bp-tfh 

We obtain after some linear algebra 

s(A, \^\)h - A. 

Setting T = 1^1 and employing the scaling z — A/r^, we obtain the boundary symbol 
s(A, r) in the form 



s(A,t) 









iiz) 



IT C 



(7.2) 



We derive this formula in the Appendix. The first and the second terms are the 
same symbols as in the system (3.1) and (3.3) in [53]. The holomorphic function 
m{z) is the same as (4.14) in [53] and £{z) is given by 



1p{z)i{z) = pifl2 



LOl{z) - 1 



^2{z) - 1 „ ^2{z) + 1 



a;i(z)(a;i(2)2 + 1) \^ ^2(2) 

^2{z) - 1 

UJ2{Z){UJ2{Z)^ + I) V ^l(^) 



P2fJ-l 



0J2{z)^ + l 

fuji{z)~ 1 



UJi{z) 



LJiizy 
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where ujj{z) = ^Ji + Pjz/Jlj and 



a;i(z)2 + l \^ W2(^) uJ2izy 



0J2iz) + l fuJi{z)-l UJi{z) + l 



Note that ujj{z) is holomorphic in the shced plane C \ (—00,——], hence the 
function ip{z) has this property in C \ {—oo,!]], with r/ — min{/ij/pj}. It is not 
difficult to see that f/; maps C+ into C+, and with V'(O) = 2(/ii + ^12), '>P{z) — >■ 
^\i ~^ ^\i ^^ 1*^1 ~^ 00 we may conclude ^(C+) C S^ for some angle 
4> < ^. By continuity of the argument function, this implies ■(/'(Szl-i-,,) C C+ 
for some 77 > 0. Therefore ^(z) cannot vanish in I]ii+^. This implies that (.[z) is 
holomorphic in this sector and in a ball B^^ (0) for some tq > 0. For the asymptotic 
of l{z) we have 

2miM2 \\fiia[p2~ VWpi) 

£(0) = 0, lim zll{z) = / ' for z e C \ ]R_. 

1^1^°° MiV/^2/p2 + M2VAil/Pl 

Thus there is a constant L — L{r, cj)) > such that 

for each cf) < ^ + rj and r < tq. Combining this estimate with the estimate for 
m[z) from Sec. 4 in [23], it is easy to conclude 

|s(A,t)|<c^(|A| + |t|) AeS5+,„ r e S,„ 

whenever 7y > is small enough. Conversely since m(0) > 0, given a small 77 > 
we find T^i € (0, vq) such that tti^z) G S^-s,, and |7ti(z)| > ^-^^^^ for all t e 5^,(0). 
This implies that there is a constant c^ > such that 

|s(A,r)| > |A| + |r|^m(0), A e S.+„ |A| < r,|r|2. 

On the other hand, choosing |A| > C\t\ we obtain 

for all A e S.|+^, r G S,,, with |A| > C\t\ and 

C > 2iaM/lpf + \co\L/lpj + |6o|/M). 
Therefore if Aq is chosen large enough this implies the lower bound 

|s(A,t)| > c^(|A| + |t|) a e S^+,„ T e E^, |A| > Aq. 

Thus this boundary symbol has the same regularity behavior for |A| > Aq as that 
for the problem in |23) . 
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The linear problem with variable coefficients but small deviations for a flat 
interface, i.e. 



dth- 



pdtu - fj.{x)Au + Vtt = pfu 
div u = gd 

lv]+c{t,x)\7^h = P^g^ 

-lfj.{x)dyvj - IfioWxWJ = gy 

-2lii{x)dywj + |7r] - aA^h = g^ 

-2lfiodyw/pl + {tt/pI = gh 

+ bit,x)-V,h/lpj = h 

m(0) = Uq 

h{0) = ho 



inlR", 




i>0 


inM", 




i>0 


onM"- 


-1 


t> 


on R"- 


-1 


i> 


on M"- 


-1 


i> 


onR"" 


-1 


t> 


onM"" 


-1 


t> 


in R", 






on R"- 


-1 





(7.3) 



can be handled by a perturbation argument in the same way as in [14. The same is 
true for bent interfaces. However, the localization argument needs some significant 
modifications, which we explain in some detail now. 



7.2. General Bounded Geometries. Here we use the method of localization. 
By assumption, dil is bounded in class C^~ and S is bounded and real analytic, 
so in particular of class C^. Therefore we cover S by A^ balls Br/2{xk) with radius 
r > and centers Xfe e S such that E n Br{xk) can be parametrized over the 
tangent space T^^S by a function 9k € C^ such that |V0fe|oo < V with ry > 
defined as in the previous section. We extend these functions 9k to all of T^^I^ 
retaining the bound on V6'fc. In this way we have created N bent half spaces Sfe 
to which the result proved in the previous subsection applies. We also suppose 
that Br{xk) C ri for each k. Set U := fl\ LijLiB,./2ixk) and Uk = Br{xk), 
k = I,.. .,N. The open set U consists of one component Uq characterized by 
dfl C Uq and m components open sets say Un+i, ■ ■ ■ , Npf+rn, which are interior to 
S, i.e., W^-^Un^i C ill. Fix a partition of unity {fk}jJo™ subject to the covering 



N+m 



{u^7^ 



'), < (yjfc < 1 and X^fli™ 'Pk 



1. Note that ipo = I 



offl, ¥>eCo-(L „ „ ___ ^^__ 

in a neighborhood of dfl. Let ipk denote cut-off functions with support in Uk such 
that i^fc = 1 on the support of ipk- 

Let 2 = {u,TT, |7r],ft,) be a solution of p.4p where we assume without loss of 
generality uq ~ ho = f = gd = 0. We then set Uk = (fikU, i^k = "ySfcTr, hk = tpkh. 
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Then ior k — 1, . . . ,N Zk — (ufc, tta;, |7rfc|, hk) satisfy the problems 

pdtUk - iJ,{x)Auk + ViTk = Fk{u, tt) in M" \ S^, 

divwfe = (V<^fc|u) in]R"\Sfe, 

■Psjufel +c(i,x)Vsfc^fe = (pkP^k9u + Guk{h) on S^, 

(7.4) 

l-^iix){Vuk + [Vuk]^) + 7rfe]j/Sfc - (T(As/ifc)i^Efc = '/'fe.g*; + Gg^ {u, h) on S^, 

- |^(.T)(Vufe + [Vufe]''')i/Efc • J^Sfe/p] + [tt/c/p] = ((5fe.g/i + Ghk{u) on S^, 

9tft-fe - {pUk ■ l^sj/|p] + b{t,x) ■ VSfc/i*;/|pl = <y5fc//l + Fh^i^) On Sfe, 

Mfe(0)=0 in]R"\Sfc, 

/ife(0)=0 onSfc, 

where 

Ffe(u, tt) ={VLpk)T^ - m(2;)[A, (^fc]u, 
Guk(h) =c{t,x){\/j:^fk)h 
Gg^{u,h) = - lp{x)(y^^(pk ^u + u® VEfc</'fc)li^Sfc - cr[Asj,,(/2fe]/ii^Sfc, 
G/ifc(u) =- I/u(a;)(Vs,(/?fe«)u + u(8)Vs,<y5fe)i^Sfc -i^Sfc/p], 
F,J/i)=(6(i,x)|(Vs,^;c)M/Ipl. 

For fc = we have the standard one-phase Stokes problem with parameters p2, 
pL2{x) on 51 with no-slip boundary condition on 90, i.e., 

P2dtUo - p2{x)Auo + Vtto = Fo{u, n) in fi \ E, 

div Mo = (\7ipo\u) in fi \ E, 

uq = on 9f2, 

uo(0) =0 in fl. 

For k ^ N + I (l = 1, . . . ,m) we have the Cauchy problem of the Stokes equation 
with parameters pi, pi{x), i.e., 

pidtUN+i - pi{x)AuN+i + Vttn+i = Fn+i{u, tt) in R", 

divMAr+; = {VipN+i\u) in R", 
UN+i{0) = inR". 

Though the right members Gu^.{h), Gg^{u,h), Gh^.{u), Fh^{h) have more time 
regularity than the corresponding data class, the terms (V</?fc)7r in Fk{u,Tr) and 
(S/(pk\u) unfortunately do not have this property. In order to remove this difficulty, 
we have to decompose the problem. Here is one major change compared to the 
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on Efe, 




inIR"\Sfe, 


(7.5) 


on Efc, 




on Efc, 


(7.6) 


on Efc, 
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construction in [14: . Consider the following problem for the functions 4)k,''pk- 

^dt4)k = dtu ■ Vipk = diY{dtUipk) 
ldvdt(t>k\ = {dtU ■ vifkl 
Aipk = div Ffe 
Id.M = iFk ■ i^l 
li>k - pdtM = 
l^k/p-dtM^O 

(|7.5p is an elliptic system for {dt4>k,'4'k) satisfying the Lopatinskii-Shapiro condi- 
tion, it is uniquely solvable and its solution satisfies V'fe, dt^ik € Lp[J\ Hp{R^\T,k)). 
Hence with (j)k{0) = we obtain by a time integration ^j. € H^{J] Hp{W^ \ S^)). 
Thanks to A(/)fe = (V(^fe|u) in M" \ Sfc we conclude that 0^ e Lp{J] H*{R" \ S^)), 
which implies 

V0fe e i/pi(J;Lp(R"))nLp(J;ij3(M»\Sfe)), Vt^^ e Lp(J; Lp(M"))). 

Defining 

Uk=Uk-V(j)k: Fk{u,TT) ^ Fk{u,TT) ~Vlpk, 

we see that c\yv Fk{u, tt) = and divufe = in M" \ S^, as well as 

{lh{u,n)\\y)^{luk\\y)^0 on S^. 
Next we define 

TTfc = TTfe - V'fc + Pdt4>k - tJ,{x)A(j)k 

and observe that 

I^fc] = Kl - [M(a^)A<^fcl = Kl - 1^(0;) (V(/), 1^)1 

l^k/pj = liTk/pj - Mx)A<l>k/pj = iTTk/pj - lp{x){W<j>k\u)/pj 

on Sfe by construction. Now z^ = (uk,T^ki iTTfclj/ife) satisfies the problem 

pdtUk ~ p{x)Auk + Vnk = Fk{u, n) 

div Ufc = 

-Psjufe] + c{t,x)\/^^hk = fkP^,9u + Guk(h) 



|-^(a;)(Vufc + [Vufc] ) -fTTfe];/- cr(As/ifc)^ = ¥'fc5fc + Gg Ju, /i) 
- |^(a::)(Vufc + [VukY)v ■ v j p\ + iTTfc/p] = Vkgn + Gh^iu) 

dthk - ipiik ■ lyj/lpl + b{t,x) ■ v^,hk/lpl - ipkh + FhM 

Mfc(0)=0 
/ifc(0) = 



inM"\ 


Sfc, 


inK"\ 


Sfc, 


on 


Sfc, 


on 


Sfc, 


on 


Sfe, 


on 


Sfe, 


in]R"\ 


Sfc, 


on 


Sfe, 




(7.7) 



INCOMPRESSIBLE TWO-PHASE FLOWS WITH PHASE TRANSITIONS 41 

where 

Gh,{u) ^Gh,{u)+2l^i{x)W^^k/pj - Mx){W(pk\u)iys, ■ v^Jp] 

FH,{h)^F,,M + {pd>^.MI{pl 

Now we know that the data in (|7.7p satisfy the assumption of Corollary 3.2 in 
[14] . Therefore the solution TTfc of (j7.7p has more time regularity, we have tt^ G 
oH^(J;Lp{n)) for each a £ (0, 1/2 - l/2p). Rewrite dLZ]) abstractly as 

LkZk = Hk + BkZ. 

Then by the same argument as in fl4J we obtain 

\\zk\\E<C{\\Hk\\F + a^\\z\\E)- 
Summing over all k, z = X]fc=i™ ^k yields 

\\z\\K<G{\\H\\r + a-'\\z\\K). 

Therefore choosing the length a of the time interval small enough, we obtain the 
a priori estimate 

\\zh<G\\H\\w. 

Since the problem under consideration is time invariant, repeating this argument 
finitely many times, we may conclude that the operator L : qE — >■ qF which maps 
solutions to their data is injective and has a closed range, i.e., L is a semi-Fredholm 
operator. 

It remains to prove surjectivity of L. For this we employ the continuation 
method for semi-Fredholm operators, which states that the Fredholm index re- 
mains constant under homotopies L{X), as long as the ranges of L{X) stay closed. 
For this purpose, we introduce a first continuation parameter a € [0, 1] by replac- 
ing the 7th equation of (|3.4p into 

dth + a(-AE)U - (1 - a) ilpu ■ lyj/lpj + b{t,x) ■ VsVM) = h on S. 

With minor modifications, the analysis in subsection 7.1 shows that the corre- 
sponding problem is well-posed for each a G [0, 1] in the case of a flat interface 
with bounds independent of a S [0, 1]. Therefore the same is true for bent inter- 
faces and then by the above estimates also for a general geometry. Thus we only 
need to consider the case a = 1. 

To prove surjectivity in this case, note that the equation for h is independent 
from those for u and tt, and it is uniquely solvable in the right regularity class 
because of maximal regularity for the Laplace-Beltrami operator. So we may set 
now h = 0. 

Next we introduce a second continuation parameter /? G [0, 1] by 

Pe,U2 = PP^Ui +Pt.9u, -(3Pt^T2{u,tt)i' = -Py,Ti{u,'k)v = Pj^g 
with T{u, n) = fi{x){Vu + [Vu]"^) - tt. 
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Again, we can prove that the a priori estimates are uniform for /? G [0, 1]. The 
remaining problem for (3^0 decouples into a one-phase Stokes problem with 
mixed Dirichlet-Neumann boundary condition in ^2, Dirichlet condition on dfl 
and so-called outflow conditions on S, and a one-phase Stokes problem with pure 
Neumann boundary condition in fli. According to [4] these are well-known to be 
solvable. This shows that we have surjectivity in the case a = 1 and /3 = 0, hence 
by the continuation method also for a = and /3 = 1. The proof of Theorem 13. II 
is now complete. 

8. Appendix 
To save space, we follow here the derivation and notation in 23^, Section 4.2. 

(a) The symbol of S^^ . To obtain the algebraic system for the symbol of [S'j^^]~^ we 
set P2 = and let k be given. Then (4.9) remains valid as well as the formulas for ai,a2. 
For Qi, Q2 we have here the equations 

— /?2 H ;Si + (piai 4- P2a-2)\i\ = \p\k 

LJ2 UJl 

— (/i2/32 + /i2Q2|C| ) + Aa2 (/ii/3i + /iiQfil^l ) — Aai = 0. 

P2 pi 

Inserting Pk from (4.9,) this system becomes 

piQi -p2a2 = 0, (8.1) 

(jiQi -I- (72Q2 = Ipjk, 
where 

Pi =PiA[ 2lfi/pj—— ;^ ^\^, ] =:piAp?, 

pi 'y(z)uji(z) ijj-i_(z) + I 

P2 = P2\[-+2lp/pj \ .""IIT! ] =■■ P2AP^- 

P2 'y(z)UJ2[z) iJ2(z) + 1 

and 

^ _piA_,_pi_ _|_ pip2'y2{z) ^ p2 0^1(2:) -1 ^^ Pi A 

1^17(2) uji{z) picji(z)(cji(2) + 1) ijJi{z)uj2{z) uji{z) + 1 ' |C|7(-z)^ 

^ P2A p2 P2Ml7l(^) , pi LJ2(2)-1 , ^ P2A 

''^ " |e|7(2) ^a;2(2) p2a;2(2)(c^2(2) + 1) a;i(2)a;2(2) ^^2(2) + 1^ "' |C|7(^)'^" 
where the scaling 2 — A/|^p is employed, and recall 



LOkiz) = a/I + Pkz/flk, Jkiz) =ljJk{z) + l/LUkiz), 7(2) = pi7l(2) +^272(2)- 

This yields the transformed interface pressures 

*i = piAai = M n„ol'„o„o l€!7(^)'fc, (8.2) 

Piy2 T"-P2yi 

P? 2 - 

712 = p2Aa2 = M n„o , „o„o lCl7(^) fc- 

Pi 92 T P2?l 

Note that p^, qq axe holomorphic in C \ (—00, — r;], r/ = min{/ife/pfc} > 0, and we have 
Pfe(O) = 1/pfc, gfe(O) = 2(pi + p,2fpk/nk, 
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and 

Pfe(oo) = l//9fc, gfe(oo) = (pi + p2)VPfcPfc- 
Therefore p\, g° are holomorphic and bounded on S,r/2+£ U ^e(O), for small e > 0. 

So we need to show that the Lopatinskii- Determinant r^{z) := Pi{z)q2{z) +P2{z)qi{z) 
has no zeros in Re 2; > 0. After some tedious algebra, expanding and collecting terms, 
we obtain the factorization r^\z) = r? (2)7-2 (2), where 

r1{z) = j{z)/[ujj{z)iJ2{z){iJi{z) + l)(a;2(2) + f)], 

and 

rUz) = i—Mz) + ^i02{z)){u;i{z) + l){u;2{z) + l){iOi{z) + l){oJ2{z) + 1) 

Pi P2 

+ 2{uji{z) - l){iJ2iz) - 1) + 2^^^(c^2(2) - f) + 2^^i^(a;i(2) - f) 

/ilP2 M2pl 

+ ^^2£1(^2(Z) + 1)('^2(2) + 1) + l^l£lj,{z)iOl{z){ujj{z) + 1) 
/ilp2 M2P1 

Obviously, rj(2) has no zeros in C\ (— 00, 0]); so we only have to look at ^2(2). One easily 
checks that each summand of r2(2) has nonnegative imaginary part, provided 2 G C is 
such that Re2,Im2 > 0, so 72(2) can only be zero if each summand is zero. But this is 
not possible, as already the first term shows. 

(b) The symbol of Sx^. To obtain the algebraic system for the symbol of [S'^^]^^ we 
set gi = and let j be given. Then (4.9) remains valid as well as the formulas for ai, 02. 
For Qi, Q2 we have here the equations 

2Af2(/32 + lCpa2) + p2Aa2 - 2^i(/3i + ICI^Qi) + PiAai = 0, 



V^;32 + \^\a2 + ^^13^ + l^lm = ll/pjj. 

UJ2 iOi 

Inserting /3fc from (4.9) this system becomes 

pioi— P2a2 = 0, (8.3) 

giai +9202 = [1/pli, 



where 



and 



pi = piA[l + 2M , ,^ , /ll^j ] =■■ P^^P"' 
■y[z)u!i{z) uji{z) + 1 

P2 = P2A[1 - 21^]--- — ^ ^^,^, ] =: P2AP2, 
7(2)a;2(2) aJ2(2) + 1 



= Pi^ ri , 1 ^1(2) -1 72(2) , ^ PiA 

^' 1^1(^1(2)7(2)^ ^a;2(2)cui(2) + l^/ii(a;i(2) + l)' '\^h(z)^^^' 

^ P2A 1 cJ2(z)-l 7i(z) 1 ^ P2A 

'^^ \^\cu2{zhizy ^ u^r{z)u;2{z) + l^ MMz) + i)^ ' \^\l{zy^'' 



This yields 



^i=piAai^[l/p] ^^ jf ,, o !Ci7(^)^J, (8.4) 

Piy2 I" P2'Ji 

712 = P2\a2 = [1/p] 0^^ o l^l^(^)^J- 

Pi 92 + pIj^i 
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As in (a) the functions p'l,q1 are holomorphic in C \ (—00,— 77], 77 — ram{ij.k / Pk} > 0, 
and we have 

pU0) = 1, g^(0) = 2(/ii+M2)V/^fe, 
and 

Pfc(oo) = 1, 9fe(oo) = (VPiMi + \/P2M2) Pk- 
Therefore pi, qi are holomorphic and bounded on S^/2+e H -Be(O), for small e > 0. 

So we need to show that the Lopatinskii-determinant r'^{z) := p1{z)q2{z) + p%{z)ql{z) 
has no zeros in Rez > 0. After another tedious algebra, expanding and collecting terms, 
we obtain as in (a) the factorization r'^{z) — r^{z)r2{z), where 

rl{z) = -i{z)/[u,{z)u^2{z){u,{z) + l){i^2{z) + 1)], 
and 

rl{z) = {uji{z) + u:2{z))(uji(z) + l){ui2{z) + 1) + 2{uj^{z) - l){u:2{z) - 1) 

+ ^[2{UJ2{Z) - 1) + {i^l{z) + 1){U2{Z) + 1)] 

/il 

+ ^[2(a;i(z) - 1) + {ojI{z) + l)(a;i(^) + 1)]. 

P2 

Obviously, rj(z) has no zeros in C\ (—00, 0]); so we only have to look at r2{z). One easily 
checks that each summand of r2{z) has nonnegative imaginary part, provided z £ C is 
such that ReZjImz > 0, so r2{z) can only be zero if each summand is zero. But this is 
not possible, as already the first term shows. 

(c) The symbol of s(A,r). We calculate the symbol s(A,r) defined by (|7.2p in the 
same way as Sec. 4-3 in (53]. It is enough to seek the solution of the problem 

^^(O) - ■ui(O) = -coi^h, 
P2{dyV2{0) + i^MO)) - pi{dyvi{0) + iCwi(O)) = 0, 

-2^29j/W2(0)+ir2(0) + 2/ii9j;™i(0)-7i-i(0) =0, (8.5) 

-2{p2/p2)dyW2{0) + M0)/P2 + 2{pi/p2)dyWi{0) - *i(0)/p2 = 0. 

Multiplying the 1st and the 2nd equations by i^, and setting Pj — i^ ■ aj for aj G C""'^, 
j = 1, 2, we have 

/32-/3i = -|?|^{(a2-«i)-coM, (8.6) 

|2/ 



y/Jl2L02l32 + yi^Ldll3i = — 151 (2|5|(/iiCli +^20:2) + {pi^/JIl/cdl)fSl + {p2^/Ti2/l^2)l32) ■ 

Combining the 3rd and the 4th equation of (|8.5p . we obtain 



a, =-2m,/?,/(^| + Mj|?|'). (8.7) 

Substituting the formula into (|8.6p and using the scaling loj — y'jrj\£^\u}j{z), 'jj = 
Pj\(.\^j{z), we solve the system for jSj 

-1 + 2/{iJi{zf + 1) 1 - 2/{uj2{zf + 1) \ fpA ^ fcoT'h 

,Mi{7i(^)-4/(a;i(z)2 + l)} ^^2{l2{z)-A/{u2{zf + l)}J\|32J \ 

where we set r = \^\ and z = A/r^. Therefore substituting the /3j and aj by (|8.7|) into 

|pw(0)] _ (piQl + P2Q2)r + {P2/tUJ2)P2 + (/Ol/rCJi)^! 

[pl " M 

we obtain (17.21). 
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